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THE  CALCULATION  OF  CONICAL  SHELLS  BY  THE  VARIATIONAL 
METHOD  OF  V.  Z.  VLASOV 

B.  A.  Konovalov,  Candidate  of  Technical  Sciences 

The  purpose  of  the  present  Investigation  is  to  obtain  equations 
for  calculating  slightly  conical  shells  of  constant  thickness  with 
allowance  for  deplanaticn  of  the  cross  sections  during  both  torsion 
and  bending  and  on  the  basis  of  these  equations  to  show  the  possibility 
of  constructing  more  accurate  solutions  for  shells  of  the  aircraft¬ 
wing  type. 

The  theoretical  basis  of  the  present  article  is  Prof.  V.  Z. 
Vlasov's  general  variational  method  of  reducing  complex  two-dimensional 
contact  problems  of  the  theory  of  plates  and  shells  to  one-dimensional 
problems.  The  entire  substance  of  the  article  is  divided  into  three 
sections,  each  of  which  is  self-contained. 

Given  in  the  first  section  is  the  derivation  of  a  system  of 
differential  equations  for  calculating  a  slightly  conical  shell  of 
constant  thickness  with  allowance  for  deformation  of  the  contour  (i.e., 
according  to. the  theory  of  moments).  These  equations  are  a  direct 
generalization  of  equations  obtained  prevlosuly  [4].  On  the  basis  of 

A  _ 
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the  results  of  the  first  section,  in  the  second  and  third  sections, 
respectively,  more  accurate  solutions  for  the  cases  of  bending  and 
torsion  of  a  simply  closed  conical  shell  are  constructed  with  the 
proviso  that  the  contour  be  lndeformable. 

Thus  the  equations  In  the  first  section  are  most  general,  and, 
what  Is  especially  Important  to  note,  on  the  basis  of  them  sufficiently 
accurate  solutions  for  a  very  wide  class  of “problems  can  be  obtained. 
The  degree  of  accuracy  of  the  solutions  will  depend  essentially  on 
the  number  of  approximating  functions  which  figure  In  the  argumentation 
and  satisfy  the  physical  meaning  of  the  problem. 

The  article  Is  devoted  mainly  to  a  study  of  the  stressed  and 
deformed  states  of  a  slightly  conical  shell  with  allowance  for  deplana- 
tion  of  the  cross  sections  during  bending  by  crosscutting  forces.  Such 
a  formulation  is  entirely  justified,  since  failure  to  take  Into  account 
the  effect  of  a  constraint  of  the  deplanation  of  the  cross  sections 
during  bending  in  prismatic  shells  can  result  in  considerable  errors 
in  the  calculation.  Professor  I.  p,  Obraztsov  (Doctor  of  Tech. 
Sciences)  studied  theoretically  and  proved  experimentally  [8,  9]  the 
fact  of  the  appearance  of  considerable  bimomental  stresses  In  prismatic 
shells  as  a  result  of  a  constraint  of  the  deplanation  of  the  cross 
sections  during  bending.  It  was  precisely  this  clrcumstanoe  that 
moved  us  to  study  the  effect  of  a  constraint  of  the  deplanation  occur¬ 
ring  during  bending  in  conical  shells. 

The  final  results  for  stresses  and  displacements  In  the  general 
form  for  any  shell  parameters  and  for  external  loads  of  sufficiently 
general  type  are  obtained  in  the  article. 

In  order  to  compare  the  solution  to  the  problem  of  torsion  of  a 
conloal  shell  .with  the  results  obtained  by  L.  I.  Balabukh  [l]  and 
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B.  P.  Tsibulya  [l4]*,  a  numerical  calculation  of  the  normal  stresses 
was  made  and  showed  qualitative  agreement  with  the  results. 

The  solutions  for  bending  and  torsion  of  straight  conical  shells 
can  be  applied  to  the  calculation  of  swept-back  conical  shells.  In 
this  case  some  of  the  boundary  conditions  must  be  written  in  a  cross 
section  along  a  slanting  edge  [8], 

There  are  no  fundamental  difficulties  involved  in  taking  the 
elasticity  of  the  embedding  into  account,  i.e.,  in  considering  the 
combined  operation  of  a  swept-back  conical  shell  and  a  center-section 
(subfuselage)  shell  [8]. 

The  proposed  method  of  calculating  for  bending  and  torsion  can 
also  be  extended  to  multiply  closed  conical  shells  [9]. 

The  author  deems  it  his  duty  to  express  his  gratitude  to  A.  N. 
Yelpat 'yevskiy,  Cand.  of  Tech.  Sci.  and  Senior  Scientific  Worker  at 
the  Institute  of  Mechanics  of  the  Academy  of  Sciences  of  the  USSR,  for 
a  number  of  valuable  suggestions  and  recommendations. 

Derivation  of  the  Equations  for  Calculating  Conical  Shells 

by  Approximating  the  Displacements  with  the  Aid  of  Power 

Functions 

Let  us  consider  a  conical  shell  of  constant  thickness  related  to 
a  system  of  coordinates  z,  S  (Fig.  l), 

where  z  is  the  longitudinal  coordinate  and  determines  the  position 
of  any  transverse  coordinate; 

S  is  the  contour  coordinate  of  a  point  in  the  plane  of  this 
cross  section  of  the  shell. 

*  Translator's  Note:  Reference  [14],  although  cited  in  the  text, 
is  not  found  in  the  list  of  references  at  the  end  of  the  article. 
Possibly  reference  [15]  was  the  one  the  author  wished  to  cite. 
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In  accordance  with  the  basic  idea  of  Prof.  V.  Z.  Vlasov's  varia¬ 
tional  method,  let  us  represent  the  longitudinal  u(z,  3)  and  the 
transverse  v(z,  S)  displacement  of  a  point  M(z,  S)  in  the  form  of 
the  following  finite  expansions: 

«(*.S)=»E  £/,  (*)  9,(5)  (*-1.2.  .  .m); 

i 

«(*.$)-£  V,(z)*»(S)  (*-1,2.  .  .»). 

* 

where  U.(z)  and  V^(z)  are  the  unknown  generalized  longitudinal  and 
x  k  transverse  f isplacement3; 

cp.(S)  and  ^k(s)  are  distribution  functions  of  the  generalized 

*■  longitudinal  and  transverse  displacements  along 
the  contour  of  the  cross  section  of  the  shell 
and  are  chosen  beforehand. 

Let  us  consider  a  system  of  3.  functions  <p1(s),  each  of  which  is 
a  power  function  of  the  order  a  with  respect  to  the  contour  coordinate 
S.  We  shall  assume  that  a  system  of  k  functions  ^(S)  consists  of 
power  functions  of  the  order  with  respect  to  the  contour  coordinate 
S.  The  approximating  functions  of  which  the  systems  cp^(s)  and  ^(S) 
are  composed  must  satisfy  she  condition  of  linear  independence  and  the 
condition  of  continuity  of  the  longitudinal  and  transverse  displacements 
at  all  points  of  the  contour. 


Pig.  1.  Overfall 
view  of  &  conical 
shell. 
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Neglecting  the  normal  stresses  acting  along  the  coordinate  lines 
S,  let  us  write  Hooke’s  law  for  the  shell  under  consideration  in  the 
following  form: 


Oil 
Oz  * 


(2) 


In  order  to  determine  the  unknown  functions  U^(z)  and  V^Cz),  let 
us  apply  the  principle  of  possible  displacements. 

Let  us  take  from  a  shell  with  the  cross  sections  z  =  const  and 
z  +  dz  =  const  an  elementary  strip  (frame)  of  width  dz.  The  isolated 
strip  (Fig.  2)  is  under  the  action  of  normal  and  tangential  stresses 
exerted  in  the  cross  sections  and  given  surfaces  by  longitudinal 
p(z,  S)  and  transverse  q(z,  S)  forces  and,  from  the  geometrical  point 
of  view,  possesses  n  longitudinal  (from  the  plane  of  the  cross  section 
of  the  shell)  and  n  transverse  (  In  the  plane  of  the  cross  section  of 
the  shell)  degrees  of  freedom.  Let  us  write  the  work  of  the  forces 
acting  on  the  Isolated  strip  during  a  possible  displacement  from  the 
plane  of  the  cross  section  uj  =  <p, (s),  when  U.(z)  =  1,  and  in  the 
plane  of  the  cross  section  v^  =  ^(S),  when  "^(z)  =  1: 


Cp  (’  +  ~y  dz)  (5)  (S )  5  as  - 

-  (j)-  ^  dz dS + (f)  P  (z,  S)  -ij  (S)  dz  dS-0  (/- 1,2...  m); 

{S)  Is ) 

(j)  (--  +  J-dz)  **  (S) 5  dS  - (p  (5)  5  dS- 

(S) 


(5) 


(S) 


S)  dtMjds)  dS 
EJ 


f  (f) ?(*.  S) ’h ( S)dzdS -> 0  (A*=  1 ,2... «). 

(S) 
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Fig.  2.  For  the  derivation 
of  the  equilibrium  equations. 
The  forces  acting  on  an  iso¬ 
lated  elementary  frame. 


Let  us  replace  the  variable  of 
Integration  S  in  Eqs.  (5)  by  s£  . 

As  a  result  of  this  substitution 
the  integrals  with  respect  to  the 
variable  contours  S  in  these  equations 
are  transformed  into  Integrals  with 
respect  to'  contours  jj  of  fixed  cross 
section  z  =  l.  Another  consequence 
of  this  substitution  is  the  possi¬ 


bility  of  comparing  the  values  of  any  function  of  the  systems  cp^(s) 
and  V'ij (S)  in  a  moving  and  a  fixed  cross  section. 

Let  us  calculate  the  integrands  in  the  function  of  the  new 
variable  s  —  . 

v 

For  each  of  the  functions  of  the  systems  q>^(S)  and  ^(s)  the 
following  equalities  hold  true: 


9<(Srh(j 


(*) 


The  derivative  of  each  function  of  the  system  9j_(s)  with  respect 
to  the  contour  is  determined  from  the  formula 


(5) 

The  linear  bending  moment  along  the  contour  of  the  shell  due  to  a 
displacement  v(z,  S)  =  |Vk  (z)  ^k(s)  is  determined  by  the  equality 

=  I  Vk(z)M„(S).  (6) 

* 

Comparing  the  displacements  in  a  moving  and  a  fixed  cross  section, 
we  obtain  a  formula  for  the  linear  bending  moment  of  the  shell  M^S) 
(Vk(z)  ■  1)  for  the  new  variable  £ 
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(7) 


M’iHt 

The  normal  and  tangential  stresses  and  their  derivatives  with 
respect  to  the  coordinate' _z  along  the  generatrices,  on  the  basis  of 
(l),  (2),  (4)  and  (5),  are  determined  from  the  formulas 


Let  us  assume  that  cpj(s)  and  Vfs)  are  power  functions  of  the 
contour  coordinate  S  with  the  exponents  X  and  q.  Since  a  and  p 
include  the  power  values  of  the  contour  coordinate  S  of  all  the  func¬ 
tions  of  which  the  systems  q?^(s)  and  s)  are  composed,  they  naturally 
pass  through  the  values  X  and  q. 

Substituting  (4)-(8)  in  Eqs.  (3),  we  obtain,  respectively,  the 
following  equalities: 

E  $  [L  x 

x  »,(s)  ?/s>  !  *  -  0  ^  u,,\  (s)  (f “ + 

'  +  (7/ '  ’] 9)  («)  + 

+  j, /; (*,  5) (s)  (ff  dzds^O  (y  =  1 ,  2  ...  m); 

(M  '  *  ' 

p*r  ’  -  2  ^(7)"  1  *  w + 

+[Z)(v;+  x 
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X<|»A  (s)lds- 


('  ^  VkMk  (7)  +  * Af*  M  dt  <t* 


{*) 


EJ 


+  <£?(*.  s)^(s)(7y+,**-0  (A- 1.2.  .  .  «)• 

w 

Neglecting  terras  higher  than  the  first  Infinitesimal  order  and 
cancelling  out  dz,  we  obtain  finally  the  following  equations: 

k 

(y= 1.2.  •  •  M);  |  (9) 

k 

(//= 1,2.  .  j 


where  7  »  £  j 


aji  **  <f>*/  (s)  ?<  (') 2  rfx; 
<*> 

bH'mj> »JW  ?',(*)  Zds’ 

U) 

Oi*" /'»/(*)<’*  (s)7'rfs; 


(*) 


<■*, (•<)?;  (s)ids-, 
io 

'**“$$*  (*)  h  (s) 2  ! 

<*\  { 
c  1  .f  .U*(«).U*  (»)«/» 

***  ~  £  V 


EJ 


(O 


(10) 


The  constant  coefficients  determined  by  formulas  (10)  are  calcu¬ 
lated  for  any  arbitrarily  chosen  fixed  cross  section. 

The  formulas  for  the  coefficients  (10)  can  also  be  extended  to 
shells  reinforced  with  longitudinal  elements.  In  this  case  the 
quadratures  (10)  must  be  understood  In  the  sense  of  Stieltjes  integrals, 
The  free  terms  in  fiqs.  (9)  are  determined  from  the  formulas 

Pj>"(j>P  (z,  S)tj(s)ds-, 

(*.  S)**  (*)<**• 


<*> 


(11) 
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Formulas  (11)  are  calculated  for  each  specific  case  of  the  distr; 
butlon  of  the  surface  loads  p(z,  S)  and  q(z,  S)  along  the  contour  S. 
For  example,  in  the  case  of  a  quadratic  variation  of  the  surface  loads 
p(z,  S)  and  q(z,  S)  along  the  contour  we  obtain 

P  (*)  “  P  s): 

?(2)  =*•*,' <7  (-•  s)  -to. 


Equations  (9)  constitute  a  nonhomogeneous  system  of  m  +  n  linear 
differential  equations  with  variable  coefficients  for  the  unknown 
generalized  displacements.  This  system  is  obtained  by  approximating 
the  displacements  along  the  contour  with  the  aid  of  an  arbitrary  number 
of  power  functions  of  any  order  of  the  contour  coordinate  S.  Therefore 
the  system  of  equations  (9)  may,  with  complete  justification,  be  cal lee 
a  general  system  of  differential  equations  for  calculating  slightly 
conical  shells  of  constant  thickness  and  arbitrary  cross  section. 
Assuming  that  each  ^bf  the  sets  of  functions  <Pi(S)  and  ^(S)  is  orthog¬ 
onal,  we  obtain 


aJi  ?,  (s)  dF  (*)  =0.  if  j#  i; 

<*) 

rMl*=j>1tk{s)h(s)dF[s)'=  0,  lf  h£k. 

(S) 


(12) 


Under  conditions  (12)  the  system  of  equations  (9)  assumes  the 

v-l.2  ...m);  I  (1J) 

r>-°  • 

(/l=al,  2  .  .  .  It). 
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For  a  system  of  1  functions  ^(s)  consisting  of  power  functions 
of  the  order  a  ■  X  ■  const  and  a  system  of  k  functi<ms  ^(S)  of  the 
order  const,  on  the  basis  of  (9)*  we  obtain 


0  -1.2.  .  .  m); 


-iHfirw+ti-yi*- o 

k 

(A  =  1 .  2  .  .  .  n). 


) 


(14) 


The  system  of  equations  (14)  for  determining  the  unknown  general¬ 
ized  displacements  is  obtained  by  approximating  the  displacements  with 
the  aid  of  a  power  function  of  any  order  a  of  the  contour  coordinate  S. 

The  system  of  equations  (l4)  is  a  generalization  of  a  system  of 
equations  [4],  where  a  =  1,  i.e.,  for  functions  that  are  linear  ^(S) 
and  constant  (of  zero  degree)  ^k(s)  with  respect  to  the  contour  coordi¬ 
nate  S.  For  shells  with  a  rigid  contour  (when  -  0)  Eqs.  (14)  are 
simplified  and  assume  the  form 


’  -£(fnv;+(ir>-° 

* 

(/  --  1, 2  .  .  .  m); 

+tf •)•>-» 

(A -1,2.  .  .  n). 


(15) 
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Equations  (15)  constitute  a  system  of  differential  equations 
with  variable  coefficients  of  Euler  type  with  respect  to  the  function. 
U±(z)  and  V^(z)  with  a  right  member  in  the  presence  of  surface  loads. 

The  unknown  functions,  which  are  the  generalized  displacements 
U^(z)  and  Vj^z)  of  the  above-mentioned  systems  of  equations  (9)  and 
(13)-(15)>  must  satisfy  the  boundary  conditions. 

The  latter  are  given  in  the  form  of  geometrical,  static,  and  mixes 
conditions,  depending  on  the-  nature  of  the  attachment  of  the  edges  of 
the  shell,  and  serve  to  determine  the  arbitrary  constants  of  integra¬ 
tion.  Let  us  show  by  examples  the  use  of  the  systems  of  equations 
that  have  been  obtained. 

A  Study  of  Cross-Sectional  Deplanatlon  During  Bending  for 

the  Case  Where  the  Deplanatlon  is  Approximated  by  Power 

Functions 

In  this  section  we  shall  seek  a  solution  for  the  stressed  and 
deformed  states  of  a  slightly  conical  shell  of  constant  thickness 
with  allowance  for  cross-sectional  deplanatlon  during  bending. 

The  deplanatlon  along  the  contour  of  the  cross  section  may  be 
given  by  one  function  or  a  system  of  functions  of  any  degree  relative 
to  the  contour  coordinate  S,  provided  that  these  functions  are  in 
keeping  with  the  physical  conditions  of  the  problem  posed. 

Let  us  study  two  problems  in  general  form. 

1.  The  deplanatlon  is  approximated  by  linear  functions. 

2.  The  deplanatlon  is  approximated  by  quadratic  functions. 

A  conical  shell  with  a  rigid  contour  (S^  *  0)  is  loaded  with  a 
concentrated  crosscutting  force  and  a  transverse  linear  load  q^zj 
given  by  a  linear  law  (Pig.  3). 
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The  concentrated  crosscutting  force  Qq,  in  keeping  with  the 
method  being  used.  Is  understood  as  a  generalized  transverse  force, 
the  effect  of  which  Is  taken  Into  account  In  the  formulation  of  the 


boundary  conditions. 

w*) 


Pig.  3.  Loading  scheme  and  geometrical 
parameters  of  a  conical  caisson. 

'/(,)  -  (*  ~  b)  |  -  e»  <  b  <  m  |; 

Cross  section  of  shell  (z  ■  I) 


The  transverse  linear  load  is  obtained  by  integrating  along 

the  contour  of  the  surface  load  over  the  possible  displacement  in 
the  plane  of  the  cross  section  of  the  shell. 


1.  The  Approximation  of  Deplanatlon  Along  the  Contour  with 
the  Aid  of  Linear  Functions 

Let  us  represent  the  longitudinal  and  transverse  displacements 
of  the  shell  in  the  form  of  the  following  finite  expansions: 

(*)*»(*)•  I 
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The  functions  cp1(s),  <p2Cs)>  and  ^i(s)  shown  ln  F1S*  4  are  chosen 
as  follows: 


*.  (S)=y(S); 


t,(5')=y(S). 


<t  i 
2 


A(S)]  +  cy(S); 


(17) 


In  formulas  (17)  the  coordinates  x(s)  and  y(s)  and  the  parameter 
dg  r&fer  to  a  moving  cross  section  z_.  The  functions  cp^(S)  and  cp2(S) 
are  linear,  while  ^1(S)  is  constant  with  respect  to  the  contour  coor¬ 
dinate  S,  and  we  can  use  the  system  of  equations  (15)  for  the  solution 
of  the  problem  formulated. 


Pig.  4.  Diagrams  of  approximating  functions 
along  the  contour  of  a  caisson. 


The  system  of  equations  (15)  will  be  simplified,  if  the  functions 
^(S)  and  <p2(s)  are  chosen  orthogonal.  The  condition  of  orthogonality 
of  the  functions  ^(s)  and  92(s)  for  a  moving  cross  section  has  the 
form 


alt 


]9i(S)9i(S)dF(S)-0. 

t) 
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Evaluating  this  condition,  we  find. the  coefficient  of  orthogonali 
zatlon  for  a  moving  cross  section 


•u. 


where  J  is  the  moment  of  inertia  of  the  cross  section  z  -  i  with 
x  respect  to  the  x-axls j 

d,,  d2,  and  F2  (cf.  Fig.  5)  are  also  calculated  in  the  cross 
section  z  *  t. 

'  • 

From  the  system  of  equations  (15)  we  obtain  for  the  given  loads 
and  displacements  (16) 

!«=.  ■£-(  '>;)- yfV'.  +  + <„V^)  -  ft 

ar^Kfr-T1 


(18) 


Equations  (18)  can  be  obtained  from  the  corresponding  system  of 
equations  given  by  Yelpat 'yevskiy  and  Konovalov  [4], 

Bearing  in  mind  the  functions  chosen  and  the  symbols  used  (cf . 
Pig*  4),  we  determine  the  coefficients  figuring  in  system  (l8)  in 
a  fixed  cross  section  z  -  1  from  the  formulas 


u) 


*11 "  K "  (*)  (IF  (s) d\  (-£•  +  +  bF^ ;  j 

(S)  '  ' 

~  -A* — (j)  f*  (s)  dF  (s)  - 
*11  ***  J*  'i?  (s)  dF  (S)  «=  2/?J; 

u) 

*i»-“  (« )  ?i  (s)  </r  (A)  =»  2c/7,; 

{•) 

w rf/r (s) ••••■  2W  +  Z7,); 


3 


rf  3 

7-7^.^“^) 


(19) 
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cn  *»<jV,  (s)  f,  (s)'c!F(s)'=2Fl-, 

<*) 

Ci3s=  C3l  *“  jjj  (s)  (s)  dF  (s)  =  2c/',; 

■<» 

'll  •*<(?❖? 

<»> 


(19) 


where  AF  is  the  area  of  the  cross  section  of  a  longitudinal  element 
of  the  shell  (longeron  belt,  stringer); 

is  the  bimoment  of  bending  inertia. 

Let  us  rewrite  the  system  of  equations  (l8)  with  the  coefficients 

(19)  as  follows: 


r  U\  +  3  U\-  7  '  ((/,  4-  cd/;  +  V',) 0; 

,»  '3  1  Vx 

//; + 3  *  u: ux -i  -f ••«/,+•  y1  tf)« os  ) 

/•'  /3  /  V  1(7,.  (7lf  t7„  /  I 

LV\  +±Ul+c(jU,2+l7U.:)-\-jV\.>.  xt  V-;«  } 

„ _ v _ f/'£.yrI _ — 


(20) 


Equations  (20)  constitute  a  nonhomogeneous  system  of  three  ordinary 
linear  differential  equations  with  variable  coefficients  of  Euler  type 
with  respect  to  U^z),  U2(z),  and  V|(z). 

Let  us  replace  the  variable  z_  by  t_  according  to  the  formula: 

/ 

Bearing  in  mind  that  U^z)  =  U^t),  U2(z)  «  U2(t),  and  Vj(z)  * 

*  VJ(t),  we  reduce  system  ('  >),  after  performing  the  necessary  opera¬ 
tions,  to  a  system  of  ordinary  differential  equations  with  constant 
coefficients 


U\ (0 + 2d/[  (0  -  NUX  (t)  -  cNU .  CO  “  'v  17 i  W  -  0; 
U\  (/)  +  2d/  ( t )  -  Id/,  (0  -  Aft/,  (/)  -  L  Vi  (0-0; 

U\  (0+ <0 +  c  \U\  (0 +  U.  (/))  +  V\ (0  +  Vx  (0- 
_ _ ssH _ (c1  —  L\ 

2 FtG{t-t)\  l)' 


(21) 


15 


whore  N,  L,  and  M  are  coefficients  determined  by  the  formulas 

Vx 


iJif  ' 

Let  us  represent  system  (21)  in  the  form  of  Table  1,  where  D 
2 

and  D  denote,  respectively,  the  first  and  second  derivatives  with 
respect  to  the  Independent  variable  t_  of  the  functions  in  the  upper 
line. 


TABLE  1 


Ui(>) 

Ut(i)  | 

! 

:  V'i  (0 

i 

flight  Month  or 

2 D-N 

""" 

*~cN 

-A' 

0 

-  L 

.  0+1 

05 +  20 -.11  —  L 

1 

0 

!  | 
c(0+l>  j  0:.  1  | 

(j  »\ 

■2flG(l-i)\‘  l) 

Let  us  introduce  into  our  study  a  new  function  f (t)  such  that: 

£/,(*)■ 


U,V). 

vko- 


~cN  —N 

D'  +  2D-M  -L 

/(O; 

D7  +  2D — N  — jV 

-L  -L  \ 

/W; 

D*-\-2D-N  -cN 

-L  D7+2D-M 

/w. 

(22) 


Then  the  first  two  equations  in  system  (21)  are  satisfied  identically, 
while  the  last  equation  gives  the  resolvent  equation  with  respect 
to. tbs  newly  introduced  function  f(t) 


(23) 


/  +5/  '-(8— P)/*  -f(4-3P)/'— 2/y'  = 

= - ! _ qJL(ct ..  t\ 

2 f,o  /- i  \  //’  ■ 


where 


P.-* 


2/Vs 


Thus  the  solution  to  system  (21)  is  equivalent  to  the  solution 
of  one  nonhomogene ous  linear  differential  equation  of  the  fifth  order 
with  constant  coefficients  (23). 

The  solution  to  Eq.  (23),  as  is  known,  consists  of  the  general 
solution  of  the  homogeneous  equation  plus  the  particular  integral. 

The  homogeneous  differential  equation  corresponding  to  (23)  has 
the  form 

/'  +  5/,V  +  (8  -P)  f"  -r  (-1  -  3P)  f"  -  2 Pf  =  0.  (  24  ) 

let  us  write  the  characteristic  equation  corresponding  to  (24) 
n-  -j-  on*  -i-  (8  -  P)  /i3  4-  (4  -  3  P)  n5  —  2P« = 0.  (  ) 

Equation  (25)  is  satisfied  for  n^  =  0.  Moreover,  Eq.  (25)  has 
the  general  roots  n^  =  -1  and  rig  =  -2  for  all  shell  parameters,  which 
can  be  verified  by  direct  substitution. 

Calculating  the  remaining  roots  of  the  characteristic  equation, 
we  obtain  finally 


where 


«i*»  — 1;  2;  n4«=  —  1  —k-  n, **--!  +  *;  nf«=0, 


We  can  now  write  the  general  solution  of  the  homogeneous  equation 


(24) 


/« (/)  + Cje-V  +  +  C4e(*-w  +  Ct. 


The  particular  integral  of  the  nonhomogeneous  eqviatlon  (23)  has 
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the  following  value j 


/<0 


00** 


M _ /. 


12^,0  (/-»)(/' -3) 


The  solution  to  Eq.  (23)  has  the  form 

/(<)-/«  (/)+/(0.  (26) 

On  the  basis  of  (22)  and  (26),  we  can  obtain  for  the  unknown 
funotlons  U^(t),  U2(t),  ard  V£(t) 

um - N  [  k'Cf-'  -t-  PC.t~*  4-  FCt + #  4- 

+  57^F<2-',<»]; 

I. 

2PjO(/-*)/»  J* 

Vl  (0  - 1 1  -f  W(/>+ 1 )+M]  C,*-'  +  WPC,e-  »  + 

-i-PlP-M)C.e-'t-",  +  P(P-M)Cie,''-i»+ 

t  A  PCj  -r 

'ijt> 


e~ 


12^,0  (/—  —  3) 


+  7 


4f,G(/-4)P 


[2(/V+A*)-W>/|. 


Returning  to  the  variable  _z,  we  may  write¬ 
up)-  -^[**C,(i)",+  />Cl(i)’’+/»C,+ 

+n^,(7r,+i?^i?(2-p'"7)]: 

j  “+,)  +  PC40.y*  '  + 


4^,6  (/ 


vo^ _ i'*A+ 1  _  Sslj  ] . 

—  *>(/>  —  3)  \ /  y  2^,0(/-*)/>  J’ 


Vi  (*)-[l  +  A/(P+  \)  +  M\Cl(jyi+NPCi(ff*+ 

+/>(p-iW)c,(i)"‘*+;,+ 

+P(P—M)Ci  (jy*~'+NPCt4 . 

4.  ftg  13(3 ~Aj)  +  A:(P_~ 3)1 ■ 

T  12/'lO(/-*)(/>-3)  w/  *  •  . 


+ 


(27) 


18. 


(28) 


Integrating  the  last  expression  in  ( 27 ),  we  obtain 
Vt  (*)«/  j[l  +N(P- 1)+ Arj  C,  In ±-NPC, 

_^(P-Af)C3(i.)-*  +  ^(P-A1)C4(i)+*  +  ArPC8(i)+,+ 

,  r  i  <?o/»[3(3-A<)  -;-*(/>  — 3))  /*\+J , 

+  ,'1'  24^,0  (/_*)(/>_  3)  W/  + 

+-,?r4^i>-  [2(A'+">  - -  •)]  (t)+'J  • 

After  the  functions  U^(z),  ^(z),  and  V^(z)  are  found,  we  proce 
to  the  determination  of  the  stresses. 


The  Determination  of  the  Normal  and  Tangential  Stresses  in 

the  Shell 

The  normal  stresses  in  the  shell  occurring,  on  the  basis  of  (2) 
during  the  displacements  (l6)  are  determined  from  the  formula 
*(z.S)—E j.V | y  C,(i.)" *+  Ct(~ }-a—  4. 


+ 


4/^,0  (/-*)(/> -3)’ ]  ^ ' 


(29) 


+• 


The  tangential  stresses  in  the  shell  occurring,  on  the  basis  of 
Hooke's  law,  during  the  displacements  (16)  are  determined  from  the 
formula 

*  (2.  S)  »<?  ( Ux  (z)  (S)  +  Ut  (2)  (S)  +  V\  (2)  *  |  (5)  | .  ( 50 ) 

A  more  accurate,  value  of  the  tangential  stresses  may  be  obtained 
from  the  differential  equation  of  equilibrium 

(5lJ 


£2{(7r"+,'s-[(T)"v'']’'<5')t+t-°- 


For  the  conical  shells  under  consideration  with  a 


1.2 


1  the 


flux  of  tangential  forces  q  ■  t6  can  be  found  by  integrating  (51): 
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(32) 


+ Tt  [(f )+  * U% ]  J ' *s  (S)  1  rfSJ +  ft  (*)• 


When  determining  the  flux  of  tangential  forces,  it  is  necessary 
to  make  a  cut  in  the  cross  section  of  the  shell  to  serve  as  the 
initial  point  in  calculating  the  integral  terms  (32).  The  flux  of 
tangential  forces  Q.q(z)  in  each  cross  section  of  the  shell  in  a  simply 
closed  oontour  is  determined  from  the  equation  of  the  moments  of  the 
external  and  internal  forces  with  respect  to  an  arbitrary  point  in  the 


cross  seotion. 

With  the  aid  of  cuts  the  cross  section  of  a  multiply  closed  shell 
should  be  transformed  into  an  open  contour,  which  makes  it  possible  to 
determine  the  integral  terms  (32).  The  fluxes  of  tangential  forces  in 
eaoh  of  the  i^  contours  qQ1(z)  are  determined  from  the  equation  of  the 
momenta  and  a  system  of  i  -1  equations  of  compatibility  of  the  defor¬ 
mations. 

The  final  expression  for  the  flux  of  tangential  forces  has  the 


forms 


«(**•$) 


0 

-ik  [(a + iiPCt^y^-Hk-i)  pc^y*-** 

wfe=*(r]f*H 


(33) 


The  transverse  displacement  (deflection),  on  the  basis  of  (i£) 
and  (aS)*  la  determined  by  the  formula 
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r  (r.  S)  ®Z  {[  1  -|-iV (P+  1 )  +  AfJ  C,  In  y  -  1  - 

_  £-(p-m)  c,  {-)"*+ j  (p-m)  ct  (±yk+NPcs(±.yi+ 


+  C.  +  ?yM3(3-Al)-!- Ar(A>-3)1  /  *  \+2 

21^,0  (/  —  i)(P  — 3)  \/J  T 


(3») 


The  Determination  of  the  Arbitrary  Constants 


For  each  particular  case  of  attachment  of  the  shell  the  arbitrary- 
constants  must  be  determined  from  the  appropriate  boundary  conditions. 

For  example,  for  the  shell  shown  In  Fig.  3  the  boundary  conditions 
have  the  form 

end  section  z  =  i  rigidly  attached: 


1)  £/.(z)  =  0,  2)t/,(z)-0,  3)  l/,  (z)  =•  0.  ' 

end  section  z  «  m  freely  deplanes: 

4)  U\  (z)  =  0;  5)  U'i  (z)  =  0, 

6)  G\cuU \  + CnUi  + rnV\\  —  —  —  Q0. 

1  I 


(35) 


Evaluating  the  boundary  conditions  (35  )>  we  obtain  a  system  of 
linear  algebraic  equations  for  the  determination  of  the  arbitrary 
constants.  This  systen  is  shown  in  Table  2.  It  should  be  noted  that 
the  system  of  equations  for  determining  the  arbitrary  constants  has 
a  very  simple  form  ahd  can  be  solved  without  any  difficulties  for  any 
parameters  of  the  shell. 

If  the  transverse  linear  load  acting  on  the  shell  is  distributed 
uniformly,  we  must  take  b  »  -«  in  formulas  (2’7)-(29),  (33)>  and  (3^) 
and  in  Table  2. 
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TABLE  2 


System  of  Equations  for  Determining  the  Arbitrary 

Constants 


*rr 

r> 

T". 

C 

O 

C‘ 

C; 

Cj  ■ 

c« 

| 

|  c  1 
1  | 
i 

C« 

j 

Rl3ht  Member 

i 

1 

i 

! 

+*J 

•  1 

1 

l 

+/» 

!  1 

‘  1 

t  j 

°  ! 

0 

_ — J9fl — (V+5J.\ 

12 FxG(t-b)P\  ^  / ; 

2 

1 

1 

-  1 

1 

! 

!  0  1 

i  i 

1 

-f  P  ; 

;  1 

i 

*r  P  1 

1 

j 

0 

j 

0 

45,0 (/- b) (P- 3) Pt  V  l  J  t  J 

3 

1 

■  0  ' 
1 

i  1 

-NP 

i 

1 

^  ( ■* 

1 

P 

: 

'  1 

1  1  ! 

-rArp’,+i 

!  40^3(3— Al)^V(P-3)l  WM2(AH-N)+A'P1 

24P,0(/- »)(P  — 3)  “  4  FiQ(t  —  b)P 

1 

4 

1 

1 

1 

+HtT  ; 

1 

0 

0 

0  \ 

B 

[/  ™\+I  «  *  /m  \+ij 

'r12P,0(/-»)K  /  /  ”  /  W  )  J 

5| 

1  j 

+  1  i 

i  1 

1 

0  1 

/  m \J 
-«+i)p(T) 

/  Jft  \4*>! 

+<*-!>/>(-,-)  ; 

! 

0 

j 

i 

0  ! 

7o (  m  \+* 

4pjO(/— t)(P — 3)  \  /  / 

! 

r> 

i 

\ 

0  i 

\ 

0 

j 

i 

0  i 

1 

0  t 

a  l 

j  On  .  7o/1  r/m\+5  7l,(m\*1] 

! 

j 

V 

iF\G  4P;  0  (/ — 6)  [\  /  )  Tvr)  J 

Approximation  of  Deplanation  Along  the  Contour  with  the 
Aid  of  Quadratic  Functions 

As  was  done  previously,  let  us  represent  the  longitudinal  and 
transverse  displacements  of  the  shell  in  the  form  of  the  following 
finite  expansions: 

u(z,  S)<=Ul(z)<tl(S)+Ut(z)<ta(S);)  (36) 

(5)-  I 

The  functions  ^(S),  V'^(s)  (cf .  Pig.  4),  and  (Fig.  5)  are 

chosen  as  follows: 


22' 


(57) 


Fig.  5.  Diagram  of 
a  quadratic  deplana- 
tion  function  and  its 
derivative  with  re¬ 
spect  to  the  contour 
of  a  caisson. 


<Pi  (>S)=y  (5); 

V.  (5)—  =F  Jf*  C^)]+cjrC5k 

(•*)=/ (S). 

The  functions  cp^(S)  and  ^(s)  remain  the 
same  as  in  Section  1.  Let  us  write  the 
condition  of  orthogonality  of  the  functions 
cpi(s)  and  cp2(S)  in  a  moving  cross  section 
a,5=^?i(5)®j  {S)  dF(S)=0, 

whence  the  coefficient  of  orthogonality  of 
the  functions  cp^(s)  and  cp2(s) 

/  6/,  ’ 

where  J  is  the  moment  of  inertia  of  the 
cross  section  z  ~  i  with  respect  to  the 
x-axis;  d1,  d2,  and  F2  (cf.  Fig.  are 
also  calculated  in  the  cross  section  z  =  i.. 

In  this  case  the  function  <P2(S)  may  be 


represented  in  the  form 


Tt(S)« 


-x’(s) 


6 J, 


y(s)  , 


(38) 


where  the  coordinates  x(s)  and  y(s)  and  the  parameter  d2  refer  to  a 
fixed  cross  section  z  -  l . 

Formula  (58)  shows  that  <p2(s)  is  a  quadratic  function  over  the 
entire  contour  of  the  shell. 

In  other  words,  the  transition  from  the  function  <P2(s)  of  a 
moving  cross  section  to  Ihe  function  <P2(s)  calculated  in  the  fixed 
cross  section  z  ■  1  can  be  made  according  to  the  formula 
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For  the  solution  of  the  problem  formulated  let  us  use  the  system 
of  equations  (9).  For  the  given  loads  and  displacements  (36)  we  have 


)  +  t“S-(vt/>)+r"S-(f,,;)- 

— _ v _ r// \+ 1 _ *_] 

G(l-b)\\  I  )  I  J  ' 

The  coefficients  a^,  bu,  b^.  -  b21,  c^,  -  Cgl,  and  ru 

figuring  In  the  system  of  equations  (39)  are  determined  from  formulas 
(19). 

For  the  coefficients  a.^2  anc^  ^22  we  ^ave  the  following  expressions 

“A,  “ 


(J) 


rfj/'t  * 
6 


A  F-i  ,1*Fi  (1F2  +  />,  +64f) 


5  37  x  ‘  367* 

=  $ <f'/ <*) <//*(«) -  • 2  d-F,-\-2c'Fl, 

%  o 

r.T» 


■} 


(*W) 


where  Is  the  blmoment  of  bending  inertia. 

After  the  coefficients  are  calculated,  let  us  rewrite  the  system 
of  equations  (39)  as  follows: 


*  2F]  u  _  *2  2c/r 


i.//  *  2^1  , , 

r/x  • 


~£/i  +  5  -  4/'-  ?-Z£LUl—*3  -hL.U  - 2££li/'-  f> 
'*  's  *  /*  771,  1  n.j,.u-  n~v'-=(y' 


1 3  77l»  •  /*  77it 


(41) 


„ _ ?r/ _ i  1 

2>,G(/- *>  [( / ;  /  [  •  I 

) 

We  have  thus  obtained  a  nonhomogeneous  system  of  three  ordinary 
linear  differential  equations  with  variable  coefficients,  each  equation 
being  of  secdnd  order. 

The  system  of  equations  (4l)  is  distinguished  by  the  fact  that, 
in  contrast  to  the  system  (20),  the  equations  figuring  in  it  are  not 
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of  Euler  type.  Consequently  we  are  unable  to  obtain  a  solution  In 
terms  of  a  resolvent  function.  The  solution  to  the  system  (^l)  will 
be  sought  by  successively  eliminating  the  functions  TJ^(z),  U2(z),  ar..'. 
V£(z)  figuring  in  it. 

Let  us  introduce  into  our  study  a  new  variable  t_  according  to 
the  formula 


-~e'. 


Performing  the  necessary  operations  on  the  unknown  functions, 
let  us  represent  the  system  of  equations  (4l)  in  the  form 
U\  (0  -r  2 U\  {() —NU,  (l)—NV\  (/)  -  (AW,  (0  =  0;  j 


-  LUX  (/)  -  LVX  (0  +  &  ((/•.  (0  +  ALT,  (0  -  MU,  (01  =  0; 

U\  (0 +£/,  (/)  +  V\  (/)  -f  V,  (/)  +  ce'lUi  (/)  +  2 £/,  (/)]  =  . 

no'- 


•2rxG(i—  b) 


(-r)- 


(te) 


In  equations  (^2)  it  is  assumed  that 


A  = 


2f,p 


/.= 


At** 


Vx 

2f  /••,/? 

T/!»  ; 
b^P 


Vi- 


The  solution  of  the  nonhomogeneous  system  (42)  consists  of  the 
general  solution  of  the  homogeneous  system  and  the  particular  Integra:' 
for  U^(t),  U2(t),  and  V^(t)  corresponding  to  the  nature  of  the  load 
under  consideration. 

Let  us  first  consider  the  homogeneous  system. 

By  eliminating  the  function  U2(t)  and  its  derivative  from  the 
first  and  third  equations  of  the  system  (42),  we  can  verify  that  the 
function  V^t)  and  the  corresponding  derivatives  vanish  simultaneously 
The  equation  obtained  will  contain  only  derivatives  of  the  function 
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U;  +  3U\  +  2U\-0. 


t 


U±(t) 


The  general  solution  fbr  U^t)  has  the  form 

£/|0,(O =C,c-»+C}e-'  +  C,.  (43) 

After  the  function  u£°^(t)  has  been  determined,  we  eliminate 
Ug(t)  and  its  derivatives  from  the  first  and  second  equations  of  the 
same  system.  The  resulting  equation  contains  the  functions  U^(t)  and 
V^(t)  and  the  corresponding  derivatives 


where 


t/r+4i/I  +  (l-iV-Af)i/I-2(3+A/+AJ)i/',+/VA'i/l- 
-  N  V ]  -  2N  V\  +  NRV\  =  0. 


R=  3  + 


3t/i» 


(44) 


On  the  basis  of  (43),  Eq.  (44)  can  be  represented  as: 

Vi + 2V\-RV\=RCxc-»  +  (l  +/?+l±^  C,e~‘  +  I<C3.  (45) 

The  general  solution  for  VJ(t)  will  consist  of  the  general  solu¬ 
tion  to  the  homogeneous  equation  corresponding  to  (45)  and  its  partic¬ 
ular  integral  and  will  have  the  form 

vr  <'»- -c-*-*-his^]c-*“'-c*+ 

40)  (t )  is  determined  from  the  first  equation  of  the  system  (42): 


cUT(t) 


e  L 


C„e--'- C<e< -*+»' '*+«)<  — Cse-<3+v 


+  1)W 

The  particular  integrals  for  U-(t),  U2(t),  and  V£(t)  in  the  system 
of  equations  (42)  will  be  sought  in  the  form 


£,W- *,*+*,/; 


(46) 
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Substituting  the  expressions  (46)  In  each  equation  of  the  system 
(42)  and  equating  coefficients  of  like  powers  ( e P^)  in  the  left  and 
right  members,  we  obtain  the  following  systems  of  equations  for  deter¬ 
mining  the  arbitrary  constants  figuring  in  (46): 

(3—  N)A,  —  cAM,  -  AM 3  ---  0;  2 B,~  cNfl,  -NB3*=0;  j 

A.+yA.-M,-.0;  (3 +  Af)B, +  £/?,-.<>,  !  (47) 

•4| -j-cA,+  ;  cB„  +  B3 •=•  +  —*2** - .  I 

4/^,0  (/—  b)  •  2/^,0  (/—»)  j 

Determining  A^,  Ag,  ...  from  the  above,  we  find  the  particular 
Integrals 


0,  (/)= _ <£* _ /V_3  *  A.  ' 

,1J  i2rlG(/-*)lc  3  /  V’ 

U..  (/)  =  --  —  (. X  JL  e  ' L.V  \  ( u  q  - 

4f,o(/-*)U  /  «~3;*  ^  (48< 

i/',(n:=r. - .*82 _ IY\-.  -^”-V'-3A/A/4.fiA,  +  3  *1 

u .  i2/,o(/-t)[r  /?-3r  3;v  / /+6— tJ-  j 

For  the  unknown  functions  of  the  nonhomogeneous  system  (42)  we 
have  the  following  formulas: 

tf,m-C1r-.  +  C^+Ci-BJte_(^- 3±,); 

—  j _ *s2£  /  1  2  b 

e  ^4FtG(/-b)\R— 3  /?  /  /’ 

v-,(0— c^-[.+ (/‘±ir]c^-.-c)+ 

+  C4C<-‘-i  1  *+T>'  -f  CgC-l1  -{- 

■f— *a2  r^-JALV'-37vA/  +  6AL+3  4  1 
Wt0  (/— S)  Lv  A-3r  ^  /  ,  +  ®~  Tj. 

Returning  to  the  variable  _z,  we  obtain 

y,(*)-c,(.fy',+c,(i)”,+c,- 
l'  (491 

w-dr^ar-wt) u,TO- 
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► 


r  (  z  v-  (2  +VW\) 

C,V77 


.+ 


V\(z) 

-c> 


.  go IJL  [  1 _ J_  _fr_  n 

r4F[0(l-t)  U-3  R  l  \l  )  J 

-c'(tr-i,+ 

+c.(f)-‘’"TO+c,(-f-) 


R  l 

4  +  A1  ] 
(tf+i)Ar 

*  x-d-M'S+i) 


■NrV 

+ 


+  6 


At  -1-3  _4_' 
R  l 


m 


The  function  V^(z)  is  determined  by  integrating  the  last  expres¬ 
sion  in  (^9) 

V,  M-<{C,  (f )■'-['  +«»^lc='"f -c*(f)+'+ 

+ys¥rc,(T)  ~7wttc,(t)  ’,”*'+c‘+ 

f  -  *)(f r + 


j  2  ilii  J.  /  *  \+ 1  ] ) 


R 


(50) 


The  Determination  of  the  Normal  and  Tangential  Stresses  in 

the  Shell 


The  normal  stresses  in  the  shell  are  determined  from  the  formula 
o(:.S)  -  /?{-  f  C,  (Ar,_-Lc!(-i-)-’-_^_[1  _ 

-3t(trih  <j')+e[-T*+vc'(fr+ 

+1^±J-C.  c.  (f) 


~<3+t'#+l) 


(51) 
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The  tangential  stresses  in  the  shell  can  be  found  according  to 
Hooke’s  law  frcm  (30). 

It  is  recommended  that  the  tangential  stresses  (or  the  Tlux  of 
tangential  forces  q  *  t6 )  be  determined  by  integrating  the  differentia 
equation  of  equilibrium  (31). 

The  flux  of  tangential  forces  in  the  case  under  consideration, 
where  =  1  and  a2  =  2,  is  determined  by  the  following  expression: 


+«fe(fr,]f’'®*«+ 


|(2  — /am- i)  v'a-r  i  c  ^  z 

(2+  //Tj-T)  /JT+T  r  (  X  (4+y/rri) 


+ 


qgbtL 


2^1  C(/ 


—  b)  u  (  /  )  ]  j  it  (S) 5  dS  J + g0  (z). 


(52) 


The  transverse  displacement  (deflection)  is  determined  from  the 
formula 


-c»(tr,+7rfrc<(  fr*'- 

^c,(f +  c.+-2i7^_  ((«_^)(f  r_ 

-^T(i"f-)(fr'+'^f(fn)M^ 


(53) 


The  Determination  of  the  Arbitrary  Constants 

If  the  end  sections  of  the  shell  are  attached  as  shown  in  Fig.  3, 
the  boundary  conditions  (35)  give  a  system  of  linear  algebraic  equations 
for  determining  the  arbitrary  constants.  This  system  is  shown  in  Table 
3. 
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For  a  uniformly  distributed  transverse  linear  load  we  assume 
in  formulas  (^9)-(53)  and  in  Table  3  that  b  ■  -  «#. 

Considering  shells  with  specific  parameters  and  load,  we  determine 
the  arbitrary  constants  from  Tables  2  and  3.  Substituting  the  values 
of  the  arbitrary  constants  and  cf  the  corresponding  coefficients  (k, 

L,  M,  N,  P,  R)  in  formulas  (29),  (33)>  (3^)  and  (5i)-(53),  we  obtain 
the  solutions  for  the  stressed  and  deformed  states  of  the  shell.  Thus 
the  problem  formulated  is  solved  in  the  general  form. 

As  can  be  seen,  the  solutions  given  in  Section  2  are  more  diffi¬ 
cult  to  obtain  than  the  solution  given  in  Section  1.  However,  it  is 
to  be  expected  that  the  stresses  represented  by  formulas  (51)  and  (52) 
will  be  closer  to  the  actual  stresses  than  those  which  can  be  deter¬ 
mined  from  formulas  (29)  and  (33).  This  circumstance  was  noted  by 
Obraztsov  [8]  for  an  analogous  problem  in  prismatic  shells. 

In  practical  calculations  it  is  entirely  permlssable  to  use  the 
solutions  given  in  Section  1. 

It  would  be  of  Interest  to  simplify  the  solutions  obtained  without 
greatly  prejudicing  their  accuracy  and  to  make  them  more  convenient 
for  practical  application. 

The  solutions  given  above  can  be  made  more  accurate  by  retaining 
a  greater  number  of  terms  in  the  expansions  (l). 

Torsion  of  a  Conical  Shell 

It  is  known  that  in  a  conical  shell  normal  stresses  appear  even 
when  the  shell  is  subjected  to  torsion  by  a  constant  moment  and  by 
free  distortion  of  the  end  sections.  The  normal  stresses  in  this  case 
are  a  consequence  of  an  Internal  constraint  caysed  by  the  conicity  of 
the  shell. 
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TABLE  5 


System  of  Equations  for  Determining  the  Arbitrary  Constants 


oi  0 

_ -  io- 

i 

0 

[ 

0 

Oo  . 

|  0  2FiG  '  4F.G(l-b)[\t  J 

1 

j 

!  *  1  j 

1  l 

!  .  -’Krf] 

Constrained  torsion  of  straight  and  swept-back  slightly  conical 
shells  was  studied  by  Balabukh  [l],  Yelpat 'yevskiy  and  Konovalov  [4], 
and  Novitskiy  [7]. 

Let  us  solve  the  problem  of  constrained  torsion  of  a  straight 
slightly  conical  shell  of  constant  thickness  with  a  rigid  contour 
on  the  basis  of  Prof'.  V.  Z.  Vlasov’s  variational  method. 

In  the  first  approximation  we  shall  assume  that  the  displacements 
of  the  shell  have  only  two  degrees  of  freedom:  one  in  the  longitudinal 
direction  (only  one  term  is  retained  in  the  formula  for  the  longitu¬ 
dinal  displacement  (l))  and  one  in  the  transverse  direction. 

( 

I 
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Fig.  6.  Loading  scheme  of  a  conical 
caisson. 

_  s»o  . 

®(I)  “  737  (*  —  1)  I  —  oo  <  X  <  m|; 

^|  *=  d\\  F*  *= 

Cross  section  of  shell  (z  =  1) 


The  shell  is  assumed 
to  be  loaded  with  a  concen¬ 
trated  torque  Hq  in  the 
end  cross  section  z  «  m 
and  a  linear  torque  jn^zj 
given  by  a  linear  law 
(Fig.  6). 

Let  us  give  the  longi¬ 
tudinal  and  transverse 
displacements  of  the  3hell 
in  the  following  form: 

« (*,  S)  -  Ux  (2)  <p,  (5);  \  /  54 ) 

*(2,5)=^,  (2)*,($)  f 


The  functions  cp^S)  and  ’^(s)  represented  in  Fig.  7  are  chosen 
as  follows: 

*,($)=/,($),  |  (55  J 

where  h(S)  Is  the  length  of  the  perpendicular  dropped  from  the  origin 
to  the  corresponding  plate  of  the  shell. 

The  function  cp^S)  is  quadratic,  while  the  function  ^(s)  is 
linear  with  respect  to  the  contour  coordinate  S. 

For  the  solution  of  the  problem  formulated  let  us  use  the  system 
of  equations  (15). 

For  the  given  loads  and  displacements  (5^)  we  obtain 


M  \(  •  \+u 
0</-X)[U)  * 


(56) 
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v^rrst? 


y 


Fig.  7*  Diagrams  of  approx¬ 
imating  functions  along  the 
contour  of  a  caisson. 


The  coefficients  figuring  in  the 
equations  of  the  system  (56)  are  calcu¬ 
lated  in  a  fixed  cross  section  z  *  1 
from  the  formulas 

(*>  dF(s)  =  y  d]d]{Fx  +  F.r\-C±F)-, 

{S) 

bn=j)  it  is  )dF(s)*=~  (d\Fl  +  (f\Ft)- 

(S) 

(s)  *1  is)  dF  (s)  =  ±  (d$F,  -  d\Ft)\ 

i*)  1 

r"  =  ?  (*)  dF  (s)  =  -L  +  d\Ft), 

(») 


where  J.  is  the  bimomert  of  torsional 
inertia. 


Let 

^ — flj|f  e 


£dW(Fx+Ft  +  teF).  | 


bt~bllG~ruG'=±-G(dlFl  +  dlF:).  \ 
b.t-cnG—\.G{d\Fx-d\F,).  \ 


(57) 


After  certain  transformations  ,  let  us  write  the  system  of  equa¬ 


tions  (56)  with  the  coefficients  (57)  in  the  form 


~  (*U\  +5 W\)—£ut~  u;  =  0; 

bt(zU\  +3  U,)+b,(zV\  +31/’,)  = 


TOq/» 

/-X 


m'-Hi-n 


(58) 


We  have  obtained  a  nonhomogeneous  system  of  two  ordinary  linear 
differential  equations  with  variable  coefficients  of  Euler  type  with 
respect  to  U^(z)  and  V^(z),  each  equation  being  of  second  order. 

We  have  already  encountered  the  solution  of  such  a  system  [system 

(20)1. 

After  replacing  the  variable  z  by  t  according  to  the  formula 


■3> 


performing  the  necessary  operations  on  the  unknown  functions,  let  us 
write  the  system  (58)  in  the  form  of  Table  4,  where  D  and  iP  denote, 
respectively,  the  first  and  second  derivative  with  respect  to  the 
indpendent  variable  _t  of  the  functions  In  the  upper  line. 

Let  us  introduce  into 


[’ABLE  4 


(0  |  V\  (/) 

1 

1  Kenber 

0  b  t 

1 

-1 

1 

!  0 

1 

+  3) 

1 

M^  +  3)  ; 

-He-H 

1 

our  study  a  new  variable 
f(t),  while  identically 
satisfying  the  first  equa¬ 
tion  of  the  system  (58) 

^1  (0=  +/; 


4/v-g/. 


(59) 


From  the  second  equa- 

tion  in  (58)  we  obtain  a  resolvent  equation  with  respect  to  f(t) 


where 


/'"  +  7/"  -}-  ( 1 2  —  A3)  /'  —  3/:1/ . 

W4h  /  ,  x 

(/— i)  V  1  c  >' 


k1-- 


jtA 

ab , 


(60) 


(61) 


Equation  (60)  is  a  nonh^mogeneous  linear  differential  equation  of 
the  third  order  with  constant  coefficients. 

The  homogeneous  differential  equation  corresponding  to  (60)  has 
the  form 

+  7 7"  4-  ( 1 2  —  k2)  /'  — 3£3/=0. 

Its  characteristic  equation  is  written  as: 

n*+  7/t* +  (12  — AJ)/t  —  3A*«=0.  (^2) 


Equation  (62)  for  any  shell  parameters  has  the  roots 
n,--3;  - 2-1/4  +  #;  «,™_2  +  i/4+p. 

Let  us  write  the  general  solution  for  the  function  f(t): 
/<«»  (0— +  C*<?-<2+,'<+F> '  +  C3f(  5  «  >  <T*  )/_ 


The  particular  integral  (60)  has  the  following  value: 

7(<)--ggf<6«  [ - 1 — c-' - 1 _ L*-*] 

(/  — X)a*t  [2(3  -4  *=)  4  +  *i  l  *  J’ 

The  solution  to  Eq.  (60)  ha3  the  form 


/(o-/«<o+7(o.  (63) 

Let  us  pass,  according  to  formulas  (59)>  from  f(t)  to  U-^t) 
and  V£(t);  we  replace  the  variable  t_  by  z_  and  obtain  expressions  for 
the  unknown  functions 


\-s+V4+F 


U,  w-  c,  (f )-+  c,  (f )- 6  Mi +  c,  ( f)- 

4.  .  f _ 1  /  *_\- 1  _  1  _x_  /■_£_ v-  2 1 

(/-X)aft,  [2{3  i  *5)W  /  4j  *  ■-•/(/  J  J; 


-(2  +  1*4+*  •> 


—  ^-C 


z  \-2*r4+* 


■(fr 


■  Sip/3  r  3a  4*  *  1  4a  +  fr|/2  X  /  *  \— 21  I 

'r(/-X)aiJl  2(3+**)W/  ^  4  +  A5  T\T)  J'  J 


(64) 


The  function  V^z)  is  found  by  Integrating  the  last  expression 
in  (64) 


v'W“,nsrc'(f)' '  + 


2ij/2 
+  : 


MH  /4-(-*J) 

2  \-l+V4+*» 


>(i) 


-0+1  i+k>) 


■(f)- 


+  + 


»l(l~/4  4-*5)  S\ 

,  TOq ft  f  3a  4- 1, ft  .  f  4a4-»|P  X  /  t  \-ni 
‘1*(f-X)at,t  2(3  +  A5)  /v  44-AJ  /  \  /  /  J|' 


(65) 
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The  Determination  of  the  Normal  and  Tangential  Stresses  in 


the  Shell 


The  normal  stresses  in  the  shell  are  determined  from  the  formula 


(i+v'i+P) 


2-/4  +  *2  n  l  *  X-3+/4+*’  L 

: - 7 — StJ  + 

The  tangential  stresses  in  the  shell  can  be  determined  from  Hooke's 

law. 

As  was  done  previously,  it  is  recommended  that  the  tangential 
stresses  or  the  flux  of  tangential  forces  q  =  t6  be  determined  by 
integrating  the  differential  equation  of  equilibrium  (21). 

In  the  case  of  the  quadratic  function  <p^(s)  (a  =  2);  after  perform¬ 
ing  the  necessary  operations,  we  have  for  q(z,  S) 

g (z,  S)=  -  £  [3C,  +(2  +  V 4+P)  |/7+S> C, X 

+(f)',,+1',~’  -(2-i/7+I=)K4+*>c,(f)-,*’OT_  (67) 

*J  J* 


The  longitudinal  and  transverse  displacements,  on  the  basis  of 
(55)>  (64),  and  (65),  are  found  from  (54).  In  practical  calculations 
for  torsion  we  are  interested  in  the  angles  of  twist  of  the  shell 
9(z)  ■  V1(z),  which  are  determined  from  (65). 


The  Determination  of  the  Arbitrary  Constants 


Evaluating  the  boundary  conditions  for  each  particular  case  of 
attachment  of  the  end  sections  of  the  shell,  we  obtain  a  system  of 
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linear  algebraic  equations  for  determining  the  arbitrary  constants. 
Let  us  consider  the  following  two  cases. 


Case  1. 

The  shell  (cf.  Fig.  6)  is  rigidly  attached  in  the  cross  section 
z  *  l  and  freely  deplanes  in  the  section  z  =  m.  The  arbitrary  constants 
are  determined  from  the  boundary  conditions 

z-l  ^ 


1.  £/,(*)« a, 

2.  VM^O. 

Zr=m 


(68) 


3.  i/',(2)= 0; 


Case  2. 

The  end  sections  of  the  shell  (cf.  Fig.  6)  z  *  i  and  z  =  m  freely 
deplane.  In  order  for  the  shell  to  be  in  equilibrium  in  the  section 
z  =  i  under  the  action  of  the  given  loads,  we  apply  the  torque 


where  -*-**  "'  ?X)(/— m)  is  a  torque  equalizing  the  linear  torque 

4  (i  *—  A/ 

JJl(2)  over  the  length  of  the  shell  L  =  l  -  m. 

The  arbitrary  constants  are  determined  from  the  boundary  condition! 

1.  z-l  U\{z)=Oi 

2.  z=m  U\(z)=0; 


3.  z=z0  ^(btu,+by\y 

4.  z—Zq.  Vt{z)~ 0. 


■H 


(69) 


i  w: 


where  Hi (f)  *=Ht  +  ~ ^  —  (z-m)-  is  the  torque  in  a  moving  cross  section 

of  the  shell. 

It  is  enough  to  fulfill  the  third  and  fourth  conditions  of  (69) 
in  any  cross  section  of  the  interval  m  £  z  £  l .  The  cross  section  in 
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which  the  fourth  condition  is  fulfilled  is  the  initial  section  when 
calculating  the  transverse  displacements. 

Let  us  write  conditions  (58)  and  (59)  in  the  form  of  Tables  5  and 

6. 


Systems  of  Equations  for  Determining  the  Arbitrary  Constants 

It  would  be  interesting  to  consider  a  shell  loaded  with  a  constant 
linear  torque  plQ. 

In  this  case  we  assume  that  X  *  -  «  in  formulas  (64 ) — (67 )  and 
in  Tables  5  and  6.  The  problem  formulated,  as  in  the  case  of  bending, 
is  solved  in  general  form.  For  the  purpose  of  comparing  the  solutions 
obtained  by  us  with  the  results  obtained  by  L.  I.  Balabukh  [l]  tod 
B.  P.  Tsibulya  [l4],  a  calculation  of  the  normal  stresses  in  a  conical 
shell  with  a  rigid  contour  was  made  for  the  following  parameter  values 
in  the  end  section  z  =  1: 

p 

d1  =  18  cm;  d2  =  60  cm;  6^  =  0.2  cm;  *  0.3  cm;  AF  *  5-5  cm  ; 
l  =  213.5  cm;  m  =  95*5  cm. 

The  length  of  the  shell  L  *  l  -  m  =  120  cm; 

The  coefficient  7  =  ^  =  2.67. 

The  conical  shell  in  the  end  sections  z  *  l  and  z  =  m  is  loaded 
with  a  torque  Hq. 

The  coefficients  (57)  and  (6l)  have  the  following  values: 
a  =  5,527  ■  104#G;  &,=9396G;  62  =  35G-lG;  *2-6G,36. 

The  roots  of  the  characteristic  equation  (62) 

n,  =  — 3;  «**=  — 10,388;  -fG,388. 

The  normal  stresses  are  determined  along  the  edge  of  the  shell 

(*($)- ) 
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from  formulas  (66).  Graphs  of  the  normal  stresses  for  the  two  cases 
of  attachment  of  the  end  sections  of  the  shell  considered  above  are 
shown  In  Fig.  8. 


TABLE  5  (Case  i) 


c 

<** 

•H 

h 

o 

*  Ct 

c’ 

C, 

C< 

4  Right  Metriber 

1 

I  !  +1 

i 

1 

i 

+1 

+1 

0 

,  ®oWj  r  -i  ,  i 

+  (/-X)a»,  L2(3+*S)  +  4  +  *s7  J 

2 

!  1 

3a  T-  |  bf  | 

+  1 

TOoX/  4«  +  »,/2 

i 

f>\  (i  —  /  4  +  &) 

+  (/— X)o»,  4+ *2 

i 

3 

+3 

{2  +  Vr^)(^yt-Vi7W 

/  x  /  m  \+i  +/4+F 

(2— /4+*J)(y) 

. 

i 

0 

i 

I 

! 

.  W*h  f  1  ("ny-1 

'  (/-X)aft,  L  2 (3  +  \  1  J  + 

,  2  X  /m  \+t1 

‘r4  +  A*  /  \  /  /  J 

4 

■ 

0 

0 

! 

i  i 

i 

0 

i 

i 

— sii(fr-Kfn 

TABLE  6  (Caae  II ) 


Condi¬ 
tion  _ 

1 

c,  j 

» 

cj  ;  c, 

\ 
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1 

+3 

2  +  /4+A2 

2  — 

.  . 

j 

0 

-1 

.2  (3  -r  *“)  + 

ft] 

2 

1 

| 

T,  /  nt  \+  1—  V4±k* 

+3  !(2  +  V7TF)(Tj 

•  j 

0 

,  w4bj  r  1  .(•)*», 

”*■(/  — L  2  (3 -h  *2)  \  l  I 

~4  +  k*  l  [i  J  J 

3 

ab , 

~V  +  V) 

do*  • 

0 

j  0 

0 

..  _/  ' 

[i(f r- 
m 

4 

3d  -r 

b2 

bt  • 

1 

+1 

TOoM  4«  +  »,/* 

2»j/J 

1 

*l(l +  /<+*>) 

+  (/— \)abi  4  +  *J 

Comment:  The  third  and  fourth  conditions  are  written  in  the  cross 
section  z  =  i . 


Curve  a  is  obtained  during  free  distortion  of  the  end  sections 
and  indicates  an  internal  constraint  inherent  in  conical  shells. 

Curve  b_  is  plotted  for  a  shell  with  a  rigid  attachment  of  the  end 
section  z  »  i .  The  difference  between  the  ordinates  of  curves  a  and 
b_  shows  the  effect  of  the  embedding  of  the  section  z  «  \  of  a  conical 
shell. 

Comparing  the  curves  in  Pig.  8  with  analogous  graphs  [l,  14],  we 
can  conclude  that  there  is  an  agreement  in  the  nature  of  the  distribu¬ 
tion  of  the  normal  stresses  over  the  length  of  conical  shells.  Such 
an  agreement  should  not  be  regarded  as  accidental. 
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As  has  been  shown  [4],  a  calculation  for  bending  of  a  conical 
shell  on  the  basis  of  Prof.  V.  Z.  Vlasov’s  variational  method  in  the 
particular  case  corresponding  to  the  hypothesis  of  plane  cross  sections 
gives  expressions  for  the  stresses  and  displacements  which  entirely 
coincide  with  the  solutions  of  the  resistance  of  materials  when  shear  is 


taken  into  account. 


Fig.  8.  The  loading  scheme  of  a  caisson 
and  its  dimensions.  Cross  section  of 
the  shell  (z  =  l).  Graphs  of  the 
distribution  of  normal  stresses  along 
the  edge  of  a  conical  caisson. 


A  more  accurate  calcu¬ 
lation  of  the  torsion  of 
a  conical  shell  by  the 
proposed  method  in  the 
second  and  subsequent 
approximations  can  be 
carried  out  by  representing 
the  displacements  (i)  in 
the  form  of  several  terms 
of  a  series,  as  was  shown 
above  in  the  solution  of 
the  bending  problem 

Thus  the  use  of  Prof. 
V.  Z.  Vlasov's  variational 
method  for  calculating 
conical  shells  makes  it 


possible  to  obtain  mare  accurate  solutions  far  the  stressed  and  deformed 
states  both  in  the  case  of  bending  and  in  the  case  of  torsion. 
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THE  EQUATIONS  OF  AXISYMMETRICAL  THREE-LAYER  SHELLS 
WITH  A  LIGHT  FILLER 

V.  F.  Karavanov 

In  a  previous  article  [1]  the  basic  equations  of  axl symmetrical 
three-layer  shells  with  a  light  filler  were  given  without  derivation 
in  a  linear  formulation  with  assumptions  based  on  neglecting  the 
longitudinally  directed  stresses  in  the  filler  and  the  bending 
rigidity  in  the  supporting  layers. 

In  the  present  article  we  shall  give  the  derivation  of  these 
equations  in  detailed  form  with  the  scheme  and  notations  of  E. 

Reissner  [2]  taken  into  account. 

The  equations  obtained  for  axisymmetrical  three -layer  shells 
are  similar  in  form  to  the  equations  of  the  theory  of  thin  homogeneous 
axisymmetrical  shells  obtained  by  E.  Meissner  [J], 

The  equations  for  axisymmetrical  conical,  cylindrical,  and 
spherical  three-layer  shells  are  obtained  from  the  general  equations 
as  particular  cases. 

In  certain  practical  cases  terms  depending  on  the  transverse, 
compressive  deformation  of  the  filler  can  be  neglected.  For  this 
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case  the  basic  equations  of  axisymraetrical  three-layer  shells  are 
given  without  taking  the  transverse  compression  of  the  filler  Into 
account. 

The  well-known  equations  for  homogeneous  thin  axlsymmetrlcal 
shells  are  obtained  as  a  particular  case  of  the  equations  of  axlsym¬ 


metrlcal  three-layer  shells. 


Basic  Assumptions 


The  wing  and  fuselage  surfaces  of  present-day- high-speed  aircraft 


must  remain  smooth  under  a  considerable  load.  There  thus  arises  the 


question  of.  the  transition  from  a  stringer  covering  to  a  covering  of 
three-layer  type,  since  the  latter  possesses  great  rigidity  and 
strength  and  yet  is  light-weight.  Great  rigidity  of  coverings  is  of 
special  value  in  the  case  of  high-speed  aircraft,  since  the  question 
of  sagging  of  the  covering  acquires  great  significance  In  connection 
with  the  considerable  increase  in  flight  speed.  Three-layer  structures 
with  a  light  filler  also  possess  good  heat-insulation  properties, 
soundproofing,  good  vibrational  characteristics,  and  other  qualities . 

A  three-layer  structure  consists  of  two  strong  thin  outer  layers 
connected  to  each  other  by  a  filler  which  ensures  the  concerted 
operation  of  the  supporting  layers.  Foam  plastics,  honeycomb  struc¬ 
tures,  corrugated  sheet  metal,  and  other  materials  may  serve  as  the 
filler. 


The  supporting  layers  are  made  of  metal,  plastic,  plywood,  and 
delta  wood.  The  calculation  of  three-layer  plates  and  shells  with 
allowance  for  deformation  of  the  filler  as  a  three-dimensional  body 
is  very  cumbersome.  Moreover,  it  is  of  little  use  in  the  case  of 
light  fillers,  the  bending  rigidity,  tensile  rigidity,  and  longitudinal 
shear  of  which  are  small  in  relation  to  the  rigidity  of  the  supporting 


layers.  In  most  cases  approximate  theories ,  based  on  various  schemes 
acceptable  In  practice,  are  used  In  calculating  such  structures. 

The  present  article  Is  based  on  a  scheme  proposed  by  E.  Reissner 
[2]  with  the  following  basic  assumptions: 

a)  the  bending  rigidity  of  the  supporting  layers  themselves  is 
not  taken  into  account  (the  nonuniformity  of  the  distribution  of 
stresses  over  the  thickness  of  the  supporting  layers  is  neglected) ; 

b)  the  filler  undergoes  only  transverse  shear  and  transverse 
compressive  deformation,  i.e.,  it  is  assumed  that  longitudinal  stresses 
and  moments  are  absorbed  exclusively  by  the  supporting  layers,  while 
the  filler  totally  absorbs  the  transverse  force  (i.e.,  that  the  moduli 
of  normal  elasticity  and  shear  of  the  filler  are  equal  to  zero  in  the 
longitudinal  directions,  but  differ  from  zero  in  the  transverse 
direction) ; 

c)  the  filler  is  regarded  as  comparatively  thick,  light,  elastic, 
and  isotropic  with  a  relatively  small  modulus  of  elasticity; 

d)  the  dissimilarity  between  the  lengths  of  the  median  surfaces 
of  the  inner  and  outer  layers  is  taken  into  account. 

The  Geometry  of  the  Shell 

We  shall  determine  the  location  of  points  on  the  median  surface 
of  an  axlaymmetrlcal  three -layer  shell  by  means  of  the  angles  0  and 
9  (Pig.  1). 

Let  0  be  the  angle  'between  the  normal  to  the  median  surface  and  the 
axis  of  the  shell; 

9  is  the  angle  between  the  two  meridional  cross  sections; 

Ri  is  the  radius  of  curvature  of  the  meridional  cross  section; 

Ra  is  the  second  principal  radius  of  curvature. 

By  means  of  the  two  pairs  of  meridional  and  normal  conical  cross 
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sections  let  us  single  out  an  element  of  an  axlsymmetrlcal  three-layer 
shell  with  the  dimensions  dSi  and  dS2  (Fig.  2)  and  apply  to  It  all 
the  Internal  and  external  stresses  (Fig.  2) . 

Let  us  direct  the  coordinate  axes  as  follows!  the  x-axls  along 
the  tangent  to  the  meridian,  the  y-axis  along  the  tangent  to  the 
parallel,  the  z-axls  along  the  external  normal  to  the  median  surface 
of  the  shell  (cf.  Pig.  3) . 

The  square  of  a  linear  element  of  the  median  surface  of  the 
supporting  layers 

dSi-dS*±+dS}x.  (1) 


Fig.  1.  Coordinate  system  of  Fig.  2.  An  element  of  an  axlsym- 
an  axlsymmetrlcal  shell.  metrical  three-layer  shell. 


The  subscript  w+"  refers  to  the  outer  supporting  layer,  while 
the  subscript  refers  to  the  Inner  supporting  layer. 


In  turn. 


(2) 


where  Al+  and  Aa+  are  the  Lame  constants  for  an  undeformed  median 
surface  of  the  supporting  layers.  , 


Fig.  3.  Stresses,  moments.  Fig.  4.  An  element  of  the 

and  external  loads  acting  outer  supporting  layer  of  a 

on  an  element  of  an  axisym-  three-layer  shell, 

metrical  three-layer  shell. 

Taking  into  account  the  dissimilarity  between  the  lengths  of 
the  median  surfaces  of  the  outer  and  inner  supporting  layers,  the 
Lam£  constants  of  a  linear  element  on  these  surfaces  are 


where 


«i± 

A 


14 


h+_t  ' 
2*1  1 


14 


h+t 


2*i 


(5) 


Here  t  is  the  thickness  of  the  outer  and  inner  supporting  layers  and 
is  taken  identical  for  both; 

h  is  the  thickness  of  the  filler  layer. 

Let  us  resolve  the  intensity  of  the  external  load  acting  in  the 
general  case  on  the'  surface  of  the  supporting  layers  into  a  force 
acting  along  the  normal  and  a  force  acting  along  the  tangent  to  the 
arc  of  the  meridian  p+  and  (cf .  Fig.  3) . 

The  third  component  of  this  intensity,  owing  to  the  condition 


¥1 


of  symmetry,  goes  to  zero 


Equilibrium  of  the  Supporting  Layers 
In  order  to  obtain  the  complete  system  of  equations  for  an 

axisymmetrlcal  three-layer  shell.  It  Is  necessary  to  consider  the 
equilibrium  of  the  supporting  layers  and  the  filler  layer  of  the 
shell  separately.  By  combining  these  equations  we  obtain  the  dif¬ 
ferential  equations  of  equilibrium  for  the  entire  composite  shell. 

The  equations  of  equilibrium  for  an  element  of  the  outer  and 
inner  supporting  layer,  respectively,  of  an  axisymmetrlcal  three- 
layer  shell  for  small  displacements  of  the  median  surface  have  the 
form 


(/V , .  alt  Rt  sin  6)'  —  Njat  ±Rcos  8  +  (p±  4- 1±)  a,  sin  6  «=  0,  (4) 


ai±#i  ai±Rt 


V±  ■+■  °±- 


(5) 


Here  Ni+  and  N2+  are,  respectively,  the  specific  normal  meridional 

t 

and  annular  stresses  of  the  outer  and  Inner  supporting  layers  of  a 
three-layer  shell  (kg/cm) ;  a+  and  r+  are,  respectively,  the  normal 
and  tangential  stresses  In  the  transverse  direction  of  the  filler 
layer  acting  on  the  outer  and  Inner  supporting  layers  of  a  three- 
layer  shell  (kg/cm2);  p+  and  q+  are,  respectively,  the  components  of 
the  intensity  of  the  external  load  acting  on  the  outer  and  Inner 
supporting  layers  of  a  three-layer  shell  (kg/cm2) ;  the  prime  indicates 
differentiation  with  respect  to  the  variable  Q. 

Figure  4  shows  an  element  of  the  outer  supporting  layer,  while 
Fig.  5  shows  an  element  of  the  Inner  supporting  layer.  All  other 
Internal  factors,,  owing  to  the  condition  of  symmetry,  go  to  zero. 

Distribution  of  the  stress  In  the  Filler  Layer 

1 

As  was  mentioned,  the  filler  is  regarded  as  light,  elastic, 
isotropic,  and  homogeneous.  Its  bending  rigidity,  tension-compression 


rigidity,  and  shear  in  the  longitudinal  direction  are  negligibly 
small.  Consequently  it  operates  in  the  transverse  direction  during 
tension-compression  and  during  shear  (Fig.  6). 

The  differential  equations  of  equilibrium  of  an  element  for  this 
stressed  state  [4]  will  be  the  followings 

[(1+i)'(,+ir)i-*  <6> 

[w'O  +-£)(!  +-£)•]+[  *-»•(»• +x)i:  (7) 

Here  a  and  r  are,  respectively,  the  normal  and  tangential  stresses 
in  the  transverse  direction  of  the  filler.  The  subscripts  £  and  Q 
indicate  differentiation  with  respect  to  the  variables  z  and  0. 

Integration  of  Eq.  (6)  gives 


(8) 


The  subscript  "0"  indicates  that  the  stresses  a  and  r  refer  to 
the  median  surface  of  the  filler  (z  ■  0) . 

Let  us  integrate  Eq.  (7),  after  having  first  substituted  expres¬ 
sion  (8)  into  this  equation;  we  obtain 


2+1 /?a  sin  6- 


j* 


i+ 


—/?,/?,  sin  9«0»«0. 


(9) 


The  resultant  of  the  transverse  shearing  stress  Q  is  obtained 
from  Eq.  (8)  in  the  following  forms 


The  Integration  extends  over  the  thickness  of  the  filler  layer 
and  also  over  half  the  thickness  of  the  surface  layers,'  as  a  result 
of  the  assumption  that  the  stresses  a+  and  r+  can  act  on  the  median 
surfaces  of  the  supporting  layers  [53.  It  follows  from  the  meaning 

u  .  J. 

of  integration  with  respect  to  +  — that  the  supporting  layers  are 
attached,  as  It  were,  to  the  filler  by  their  median  surfaces.* 

From  Eqs.  (8)  and  (10)  It  follows  that 

(11) 

k+t 

u  I  ^ 

This  equation  was  obtained  by  discarding  the  terms  — ^ —  In  com¬ 
parison  with  unity. 

From  Eq.  (11)  It  follows  that  the  transverse  tangential  stresses 
are  uniformly  distributed  over  the  thickness  of  the  filler. 

From  Eqs.  (8),  (10),  and  (11)  the  following  dependences  can  be 
obtained: 

2"  (fli+aa+t+  +  )*■  Q; 

(12) 

o 

ai+at+~+  “  —  *£-  •  ( 13) 

When  h  +  t/2R  «  1,  using  Eqs.  (9)#  (10),  and  (11),  we  obtain 

/?,/?,  sin  6  (a1+a,+o+  -  -—(/?,  tin  6Q/;  ( 1^) 

ai+°a+®+ + .  (15) 


*  Although  this  Introduces  a  certain  Inaccuracy,  it  Is  so  small 
that  It  has  no  significant  effect  cm  practical  calculations  In  the 
case  of  comparatively  thick  fillers  forwhlcht/fc  <  0.1.  The  error 
of  this  assumption  is  all  the  greater,  the  greater  the  ratio  between 
the  thicknesses  of  the  layers  t/h.  S  •  '■ 


Equilibrium  of  a  Composite  Shell 


In  order  to  obtain  the  equations  of  equilibrium  of  a  composite 
three-layer  shell,  let  us  Introduce  the  resulting  expressions  for  the 
external  and  Internal  stresses  of  a  three-layer  shell  [2], 

The  specific  normal  stresses 

JV,-a,+;Vi++a,_/V1_;  M,=a1+/V,+  +  a,JVa_.  (l6) 

The  specific  bending  moments 

Af|-^(«I+/Va+-eIJVl_).  (17) 

The  components  of  the  external  force 

P,mai+at+P+  +  <} = (18) 


The  Intensity  of  the  moment  formed  by  the  action  of  tangential  surface 
forces 


m 


h±t 

2 


(ai+at+P+  ai-at t-P-), 


(19) 


and  also 


g 


'■j  (fll+a*+?+  )• 


(20) 


To  the  above  relations  let  us  add  the 
previously. 

The  specific  transverse  force 


expressions  obtained 


*  j  (ai+a*+t+  +  |  (12) 

The  normal  transverse  stress  In  the  median  surface  of  the  filler 


j  (flj+a*+®+ + (15) 

Using  Eqs.  (12) -(20),  let  ud  add  and  subtract  Eqs.  (4)  and  (5). 
After  certain  transformations  we  obtain  the  equations  of  equilibrium 


of  composite  axisymmetrical  shells 


(NtR3  sin  8)'  -  Ay?,  cos  8  +  (/>  + -£-) /?,/?, sin  8.0; 
(M,Rt  sin  8)'  -  Afjff,  cos  8  +  («  -  Q)  *,/?,  ,m  8-  0; 


_ L_ 

A,  RiPjtla  9 


(Q^smey— o; 


(21) 

(22) 

(23) 

(24) 


The  latter  equation  becomes  an  identity,  when  the  bending  moments 
are  expressed  in  terms  of  the  stresses. 

From  Eqs.  (9)  and  (10)  it  follows  that 


RA  sin  6o = /?,/?,  sin  e3o--i-  (Q/?,  Sin  6)'. 


(25) 


According  to  Eqs.  (23)  and  (25) >  we  obtain  the  law  of  change  of 
the  transverse  normal  stress  over  the  thickness  of  the  filler 


■  z  t  N\  ,  n3  \ 
‘  *+<(*!  +  f,  V 


(26) 


Hence  it  can  be  seen  that  the  transverse  normal  stress  consists  • 
of  two  terms,  one  of  which  is  constant,  while  the  other  varies  linearly 
over  the  thickness  of  the  filler  layer. 

The  Potential  Energy  of  Deformation  of  the  Shell 
The  potential  energy  of  deformation  of  an  axisymmetrical  three- 
layer  shell  is  the  sum  of  the  potential  energies  of  the  deformations 
of  the  supporting  layers  and  the  filler  layer 

n-n.+n,.  (27) 

The  supporting  layers  are  characterized  by  the  following  elastic 

constants:  the  modulus  of  normal  elasticity  the  Poisson  ooef- 

Fu 

flclsnt  |x  •  and  the  shear  modulus  0g  ■  • 


The  elastic  constants  of  the  filler  are  as  follows:  the  modulus 

of  transverse  elasticity  E  :  the  Poisson  coefficient  u„  Is  taken 

c  c 

equal  to  zero  on  the  basis  of  the  assumption  made  concerning  the 

operation  of  the  filler  layer;  then  the  shear  modulus  of  the  filler 
Ec 

°c  "  T“  • 

The  supporting  layers  of  the  shell  being  considered  undergo  only 
tension-compression  deformation;  therefore  the  potential  energy  of 
deformation  of  the  supporting  layers 

— 8  d8  dr  +  — 

— 2pNl_N1_)  sin  0  d8  dr. 


From  Eqs.  (l6)  and  (17)  we  obtain 
2aa+A^»^  +  -j~Af.;  2 

2a<-A/5-==A^-777Af*- 


(29) 


Using  the  relations  obtained  and  the  condition  h  +  t/2R  «  1, 
expression  (28)  will  assume  the  following  form: 


2  2|»AT,/V,)  + 

/?,/?,  cm  e  </•  ^ , 


(50) 


where  C*  is  the  tension-compression  rigidity  of  a  three-layer  shell; 

D*  is  the  bending  rigidity  of  a  three-layer  shell 
* 


C*~2EJti 


(51) 


On  the  basis  of  the  assumed  distribution  of  the  stress  In  the 
filler  layer  (cf.  Tig.  6),  the  potential  energy  of  deformation  of 
this  layer  will  be  as  follows: 


n‘  "T  fjj(£ + itX1  +j^)(1 

F 

As  previously,  we  neglect  the  terms  z/R  in  comparison  with  unity, 
while  the  values  of  a  and  r  are  taken  from  Eqs,  ( 11)  and  (26), 

Then,  Integrating  with  respect  to  z,  we  obtain 


+7M^r)/?*/?jSin9rfe^- 


(32) 


On  the  basis  of  (27),  (30),  and  (32),  the  potential  energy  of 
deformation  of  an  axisymmetrlcal  three-layer  shell 


n“T  if  F(/V*  + Wf-WVy+i  (Mj-fMj— 2|»Af,Af,)  + 

+ Hr  H+s(^+^~ »)’] +^1  •** 


(33) 


Relationships  Between  the  Stresses,  Moments,  and  Displacements 
In  the  Case  of  a  Composite  Shell 

According  to  Castlgllano's  principle,  the  stressed  state  actually 
arising  In  the  shell  differs  from  all  statically  possible  states  In 
that  the  potential  energy  of  deformation  of  the  shell  H,  determined 
by  Eq.  (33)*  assumes  a  minimum  value. 

Thus  It  Is  necessary  to  find  the  minimum  of  a  functional  H  depend¬ 
ing  on  six  unknown  functions  Ni,  N2,  Mi,  Ms*  Q*  and  c0*  which  satisfy 
the  four  equilibrium  equations  (21) -(24).  The  problem  of  finding  this 
arbitrary  extremum  of  the  functional  II  in  the  calculus  of  varia¬ 
tions  Is  usually  replaced  by  the  problem  of  finding  the  absolute 
extremum  of  a  certain  other  functional  r. 

In  order  to  set  up  the  functional  r,  let  us  use  the  method  of 
undetermined  multipliers,  which  preserves  the  complete  equivalence 


of  the  variables  [6] 


Mj,  M j,  Q,  ot)d9d<f. 


m 


In  turn 


where  P  is  an  integrand  expression  for  the  potential  energy  II  and  is 
determined  by  expression  (33);  (0,q>)  (i  *  1,  2,  3#  *0  are 

unknown  nonvarying  Lagrange  multipliers; 
v1  are  the  corresponding  equilibrium  equations  determined  by  the 
expressions  (21) -(24). 

It  can  be  shown  that  each  of  the  four  nonvarying  Lagrange  multi¬ 
pliers  is  nothing  other  than  a  displacement  in  the  corresponding 
directions  =  u,  v2  =  p,  v3  =  w,  v*  =  k, 

where  u  and  w  are  linear  displacements  of  the  median  surface  of  a 

three-layer  shell  in  the  direction  of  the  x-  and  z-axes; 
3  is  the  angle  of  rotation  of  the  normal  to  the  median 
surface  of  the  shell  in  the  direction  of  the  x-axis; 
k  is  a  quantity  proportional  to  the  transverse  compressive 
deformation. 

Thus  Eq.  (21)  Is  multiplied  by  u,  Eq.  (22)  by  p,  Eq.  (23)  by  w, 
and  Eq.  (24)  by  k. 

Then  the  functional  <&  will  be  written  thus: 

0  “  2  + -5r (^ + AfI-2MV>f*>+ 


+#  ((Ay?,  sin  ey  -  Ay?,  cos  e+(p+ sin  e] + 
+ P  KM  sin  *y  -  MiR,  cos  9 + («  -  Q)  RtRt  tin  6] + 
+  w[(Q/?Jsin6)'-(^l.+-g.-^/?,/?!l8in«]+ 


The  extremum  of  this  functional  occurs  for  the  Euler  condition 


(56) 


where  is  the  partial  variation  of  $  with  respect  to  Ntj 

is  the  partial  variation  of  ®  with  respect  to  etc. 

“3^ 

Evaluating  Eqs.  (36) ,  after  elementary  transformations,  we  obtain 
the  following  expressions: 


r,  1  (»aoc* 

i  m,  r  (a  +  oc* 

]  (h  +  t)q  . 

l*  1  >2 £<*’  , 

J  C*  12ECR,R3 

1  *. 

r  12£^,  ' 

(37) 

",  L  _  _(*  +  <)  C*1 

Hoot  6  + w  ,  (*  +  *)?. 

(38) 

l1+  12£,*i  J 

€•  r  12£«#,/?,J 

1  *» 

12 £«#,  * 

B' _ k 

D •  R\  (A-rO*i*l*»*  ’ 


- - - -  —  cetfl 


<*  +  * 
w'  —  a 


<*  +  /)0« 


*1 


Ec 


RtRxHn t(k  +  t)  * 


(39) 

(40) 

(41) 

(42) 


Let  us  reduce  the  system  of  equations  (37)  **(42)  to  a  more  con¬ 
venient  form. 

Let  us  introduce  the  variable  coefficients  Xi,  Xa,  and  Xs 


1 

(*+0/ 

£« 

2 

; 

£« 

1 

(*  +  /)/ 

£. 

7 

*« 

1 

(*  +  «< 

A. 

I 

(45) 


Then  Eqs.  (37)  and  (38)  will  assume  the  form 

(1  +j^)k-(f-}»,)Ar1-<>[JWtt±J),+^.f]. 


Let  us  eliminate  k  from  Eqs.  (39)  and  (4o)  by  using  relations 
(24)  and  (42);  we  obtain 

d+yjlf.-fr-y*.  D*(i.+1tr) 

In  transforming  Eqs.  (39)  and  (40)  it  is  assumed  that  Xi  ■ 

“  (h  +'  t)E  R'^ —  on  ^e  basis  of  the  expressions  for  D*  and  Xx  from 
(31)  and  (43). 


As  a  result  we  obtain  a  system  of  equations  for  an  axisymmetrical 
three-layer  shell:  three  equilibrium  equations  (21) -(23)  and  five 
relationships  between  the  stresses  and  the  displacements  (37) -(4l). 
Prom  these  eight  equations  eight  unknowns  are  determined:  three 
stresses  Nj,  N2,  and  Q,  two  moments  Mi  and  M8,  and  three  displacements 
U,  w,  and  p. 

The  normal  transverse  stress  in  the  median  surface  of  the  filler 
a0  can  be  determined  from  Eq.  (24) .  The  action  of  the  transverse 
shear  of  the  filler  figures  only  in  Eq.  (4l) ,  while  the  action  of  the 
transverse  compressive  deformation  is  reflected  in  Eqs.  (37) -(40). 

If  in  Eqs.  (21) -(23)  and  (37) -(4l)  it  is  assumed  that 

00,  C*— £A 

and  m  «  0>  we  .obtain  a  system  of  equations  for  thin  homogeneous 
axisymmetrical  shells.  Moreover ,  these  equations  completely  coincide 

with  the  well-known  equations  [7]. 

»  % 

Reduction  of  the  System  of  Equations  to  Two  Simultaneous 

Equations 


Let  us  reduce  the  system  of  equations  obtained  (21) -(23)  end 
(37) -(41)  to  a  system  of  two  second-order  differential  equations. 

From  Eqs.  (40)-(4l)  let  us  determine  the  specific  bending  moments 


Here 


8-(l+X,)(l+k)-(|i-X«)*. 


(*5) 


Let  us  introduce  Meissner*  s'  function  V  »  QRa.  (^6) 

% 

Let  us  introduce  the  values  of  Mi  and  M2  from  (44)  into  Eq.  (22) . 

As  a  result,  after  certain  transformations,  we  obtain  a  deforma¬ 
tion-compatibility  equation  containing  the  variables  p  and  V 


•^-ootle+(-^-7l+x;7j)oort  6+t,Jp-8Ml 


(*7) 


Here 


Ti1 


i — x. 


Tj 


«L±2l 

i+x, 


T»* 


1 

1+*.’ 


(48) 


In  order  to  obtain  the  second  equation,  let  us  consider  the  con¬ 
dition  of  equilibrium  of  the  part  of  the  shell  cut  off  by  a  parallel 
circle  of  radius  r  «  R2  sin  Q  (Fig.  7),  i.e.,  let  us  project  all  the 
forces  acting  on  the  part  of  the  shell  and  on  the  axis  of  symmetry  of 
the  shell.  Then 

S00i=»2it/?j$ln6(A^,sinfi—Qco«8)  —  2*C-f 
+  J  2*8,  stn  8  (/>  sin  9  -  q  cos  8)  /?,</« -0. 


This  equation  can  be  conveniently  rewritten  as: 

F  (8)  ■■  ( sin  8 — Q  cos  0)  Rt  sin  8  m 
C+  J  *!n  8  (q  co*  8 — p  tin  6)  /?,  . 


m 


<4 


Here  F(0)  is  a  function  of  the  load  and  the  geometrical  dimensions 
of  the  shell  and  is  the  axial  component  of  the  internal 
forces  acting  on  one  radian  of  the  contour  of  the  cross 


section; 

2irC  ■  Po  is  the  axial  component  of  the  concentrated  forces  (not 
shown  in  Pig.  J) . 

Let  us  express  Ni  and  Nz  in  terms  of  V  and  F(0) . 

Then  from  Eqs.  (46)  and  (49)  we  obtain 


1  O  .1.9 


~oot 


From  Eqs.  (23),  (46),  and  (50)  we  obtain 


(50) 


N, 


V' 


F(t) 


r  __  _ • 

1  Ri 


(51) 


Fig.  5.  An  element  of  the  Fig.  6.  Model  of  the  stressed 

inner  supporting  layer  of  a  state  of  the  filler, 

three-layer  shell. 

% 

Let  us  introduce  the  second  load  function 


N(9) 


_  F(t) 
tin*! 


+■  <}RiR,- 


Equation  (51)  will  assume  the  form 


(52) 


Equations  (21)  and  (23)  are  identical  to  zero,  if  we  substitute 
in  them  the  values  of  Nx  and  N2  from  (50)  and  (53) . 

Let  us  eliminate  w  from  Eqs.  (37)  and  (38).  We  obtain 


+  (l  —  «oot6). 


(54) 


Let  us  substitute  in  Eq.  (54)  the  value  of  u  from  (4l)  with  (46) 


taken  into  account.  We  obtain 

{". [(» +JX')R'  +(“  -T**)*]- 

-Ni  [(p -jX,)/?,+(l  +7^)^]}+ W — 

- ™L  » 

(*+/)Oe  *,  • 


(55) 


Let  us  represent  Eq.  (38)  in  the  following  form: 

Differentiating  once  both  sides  of  this  equation  and  substituting 
in  the  equation  obtained  the  value  of  u  from  (41)  with  (46)  takep 
into  account,  we  obtain 

U-w'aoi  =  +7Xl)^*" (p— + 

-“•!»)  (56) 

_*+L*a8a'+fB- _ - _ 

12ft  ¥  |/  (A  «j- 1)  OefljJ  *1»  28 

Equating  the  left-hand  sides  of  Eqs.  (55}  and  (56)  and  substi- 

i 

tuting  in  the  equation  obtained  the  values  of  Nt  and  Na  from  (50)  and 
(53)  and  their  derivatives  N(  and  Na,  after  a  number  of  transformations 
we  obtain  the  second  equation  (the  equilibrium  equation)  * 


+?x^}v’-[(1+7l0i“‘’,+ 


*1 

c* 


(A +  00*  Ri 

otft0]  v+^c*P-oa(8). 


(57) 


Here 


where 


. . Rt  *|n*  9tan«  *  • 

D> W-(t  -l-j)-.)/?,  +(«t-JX3)/?*+(!‘-Jx»)«i1*9 


—  *j  *-3  ft}  tar*  09 

*,(»>  “(«*-}  *»)#.+(«  +{>-,)/?,+(!  +ix,)AUnO  +  y 


t«rf* 


(58) 


Thus  Eqs.  (47)  and  (57)  are  the  basic  differential  equations  of 
axisymmetrical  three-layer  shells  in  the  case  of  small  displacements. 

Assuming  in  Eqs.  (47)  and  (57)  that 

Gc*-£e— oo,  ).,-»X1-.X1=0, 
ri  — H.  C*”’fA*nd,,l“‘0’ 

we  obtain  the  general  equations  of  E.  Meissner's  theory  of  thin 
elastic  homogeneous  axisymmetrical  shells  [3]. 

'  After  solving  Eqs.  (47)  and  (57)  for  the  given  axisymmetrical 
three-layer  shell  and  the  appropriate  boundary  conditions,  let  us 
determine  the  specific  normal  stresses  Na  and  Na  from  formulas  (50) 
and  (53),  the  specific  bending  moments  Mi  and  Ma  from  formulas  (44), 
and  the  transverse  force  Q  from  expression  (46) .  The  specific  normal 
stresses  in  the  outer  and  inner  supporting  layers  are  determined  on 
the  basis  of  expression  (29) .  Since  the  thicknesses  of  the  outer  and 
inner  layer  are  identical,  the  normal  stresses  in  these  layers  will  be 
equal  to  the  specific  normal  stresses  divided  by  the  thickness  t. 


(f! 


The  transverse  tangential  stress  In  the  filler  Is  calculated 
from  formula  (11),  while  the  transverse  normal  stress  Is  determined 
from  Eqs.  (24)  and  (26).  '  " 


Particular  Cases 

From  the  basic  differential  equations  for  axlsymmetrlcal  three- 
layer  shells  (47)  and  (57)  we  obtain  the  equations  for  certain  types 
of  shells:  conical,  cylindrical,  and  spherical. 

Conical  Shell 

In  the  conical  shell  the  angle  0  Is  constant,  while  the  radius 

of  curvature  Rx  ■  «o.  Let  us  Introduce  a  new  Independent  variable  £, 

the  distance  from  the  apex  of  the  cone  along  Its  generatrix  (Fig.  8) . 

We  have  dsx  =  Rxd0. 

Let  then  ^ 

ds  dt  ds  dp  dfl 

Equations  (47)  and  (57)  will  assume  the  form: 


+ [(*§■) + '*Hoot  •+***»-■£■  *«)]  p)  - 

—  [l* «xrt  3 6  4- +  -hXj^yKrt  e+T,jx 

XP-T»2*&-V-0,(S); 


(• +MJwr+K«  •+*(£)! 

*]v+- W-0,*(s). 


(59) 


(60) 


Here 


~(1'  t+iW]  £— (•jjf'+T.)#.  -£--aTt-^7L»i. 


(61) 


where 


£■(*)— (l  +  jxj)^i  +(j*  — 

1  +ix,)/?t+(i  +Ix,)/?l^ta8+  (62) 


For  a  conical  shell  it  is  necessary  to  assume  that 
/?!  =  CO,  /?,«$«*  9,  Xj«X3««0,  >.5=i23*£t 
2 


.-■L  «♦')<& ,  e-i+w.  Tl - ' 

Cc  , ,  »w< 

, 

Tt"'T»*r_ 


(63) 


I+- 


Dividing  both  sides  of  Eqs.  (59)  and  (60)  by  Rx  and  substituting 
expressions  (63)  in  these  equations,  after  certain  transformations, 
we  obtain  the  basic  equations  for  conical  three-layer  shells , 


1(1— li*)s4+(2— nI)«»t«?6s,+  0]*sp  + 

+ K1  —  t1*)  s4  +  (2 + 11)  ®ss + »W  0]  p — [(1  —  f»J)  4*-j- 

+  (1  —  2j»)  «>W  6]  «P  —  (( 1  —  t»J)J  S4  -f  2  (1  —  jxs)  «*an*  9s*  + 

+  « w  ej  —t**  9  —  0*  (5); 


(») 


+  C*Ptm0»«O?(s). 


(65) 


Here 


O’,  (s)  -  [(1  -  v?)  s*+ •  W  ® -  2H<W  9)^2-  sg- 
-10-1 **)  a*4-  •*»*  ®I -T-  *V-*[0- 1**)*  «4  + 

Cc 

-|-2(1— iis)<rt«rf9s*+«J*w»40]~  , 


(66) 


If  it  is  assumed  in  Eqs.  (64),  (65),  and  (66)  that  E_  -  CL  - 

c  c 

Eh® 

»  »,  a>  m  m  ■  0,  D*  «  ~Y*r  >  and  C*  ■  Eh,  we  obtain  the  well-known 
equations  for  thin  homogeneous  conical  shells  [8], 


Let  us  divide 
equations 


Cylindrical  Shell 

Eqs*  (59)  and  (60)  by  Rx  and  assume  that  in  these 


<*  +  <)<  E, 

2a*  Et  ’  '*  1+i,  ’ 


\+h' 


where  x  is  the  distance  from  the  left  edge  of  the  cylindrical  shell. 
We  have 


oK 


O+*j)0* 


V—-- 


H-V-V 


(l  +  ^j)0* 


ma. 


(,+T*’) 


aV. 


(*  +  0«O 


-0  +  C*p=-a*(^  +  ^ 


(67) 

(68) 


If  it  is  assumed  in  Eqs.  (67)  and  (68)  that  E  »  Q  ■  •*,  D*  » 

c  c 

2^3 

*  C*  *  2h,  ana  \8  »  m  ■  0,  we  obtain  the  equations  for  thin 

homogeneous  cylindrical  shells  [8]. 

Spherical  Shell 

Let  us  assume  that  in  Eqs.  (47)  and  (57) 

const,  a=. 

*  2  a*  £t  ’ 

7,-1. 

Then  we  obtain  the  equations  for  spherical  three-layer  shells 


r+~t  v- 


b*  * 


(l+^-tK-X)*  a* 
1+X  D> 


(69) 


V"+oot  W-  «rf?«V+  3 


[  (,+ix)<*+,)°«  J 

r^f-TTT  —’I*’-/1  +r£-V’-- 
3X  [“(l  +  71)*  J  V  ,+W 


(70) 


If  it  is  assumed  in  Eqs.  (69)  and  (70)  that  E„-  G„  ■  «*, 

.  c  c 

>•— fft— 0,  P*=-~-».na  C*— £A, 

we  obtain  the  well-known  equations  for  thin  homogeneous  spherical 
shells  [8].  Equations  (67) -(70)  coincide  with  the  corresponding 
equations  of  E.  Reissner  [2], 

The  Equations  for  Axlsymmetrical  Three-layer  Shells  Without 
Allowance  for  Transverse  CongreVslve  Deformation  of  the  filler 


In  solving  the  problems  of  transverse  bending  and  over-all  loss 
of  stability  of  three-layer  plates  and  shells,  the  effect  of  trans¬ 
verse  compressive  deformation  of  the  filler  is  generally  neglected. 
From  Eqs.  (47)  and  (57)  it  can  be  seen  that  the  effect  of  transverse 
compression  of  the  filler  may  be  neglected,  if 

2  £t  ^ 

Abiding  by  this  criterion,  we  shall  neglect  terms  depending  on 
the  transverse  compressive  deformation  of  the  filler. 

Then,  assuming  that  in  the  general  equations  for  axlsymmetrical 
three-layer  shells  (47)  and  (57)  E„  -  Xi  -  X2  -  X3  ■  0, 


we .  obtain 


Ti-P.  7i  — V- 1«*  a— 1—  p*. 


1  V+[A.„  •+(-*-)']  V- 

-irrior  V+Ric,f  ■"<*>' 


(71) 


(72) 


Here  D 


B# 

m5 


Is  the  cylindrical  rigidity  of  a  three-layer  shell. 


In  operator  form  Eqs.  (71)  and  (72)  will  be 


(73) 


Lm-rt-^V - Mim; 


(71*) 

(72*) 


where  the  operator 


Let  us  obtain  the  equations  for  the  particular  cases. 

Conical  Shell 


Assuming  that  in  Eqs.  (64),  (65),  and  (66)  Ec  -  * 

and  0  m  0,  we 

obtain 

(75) 

,  V-  +  v — [1  +  -  ~  w  fl]y+opf>6-/e  (<), 

s 

(76) 

where 

"V*  *•»  •»-***•• 

(77) 

In  operator  form  Eqs.  (75)  and  (76)  will  be  written; 


L'  (V)_"£77 Jo?*1' 9  -7+ ^  9P  <*). 

where  the  operator 

M  )-*<  )•+(*)•— ! y-.  (78) 

The  solution  pf  Eqs.  (64)  and  (65)  for  a  conical  three-layer 
shell  with  transverse  couqpressive  deformation  .of  the  filler  taken 
into  account  gives  rise  to  serious  mathematical  difficulties.  On 
the  other  hand,  Eqs.  (75)  and  (76)  for  a  conical  three-layer  shell 
without  allowance  for  transverse  compressive  deformation  of  the  filler 
have  a  simpler  mathematical  appearance. 


Cylindrical  Shell 

If  it  is  assumed  in  Eqs.  (67)  and  (68)  that  E„  -  «*  and  Xa 

C 


0, 


we  obtain  the  equations  for  cylindrical  three-layer  shells  without 
allowance  for  transverse  compression  of  the  filler 


at _ 

P"  d  o' 


,v“-ii+Wc*+c*?-  -(*•+ ■?•)•'• 


(79) 

(80) 


Spherical  Shell 

In  the  case  of  a  spherical  shell  we  assume  that  in  the  general 
equations  (71)  and  (72)  Rx  «  r2  «.  a  (where  a  is  the  radius  of  the 
sphere) .  We  obtain 


w {V) ~kk+ko; *] v+ -«v-o +i»)m 


(81) 

(82) 


0? 


JTD-TT-62-1652/1+2+4 


Here 


m  )■(  r+(  )'«»*•-(  )«•**•, 

n  (j  +  »o* 

Boundary  Conditions 


(83) 


In  order  to  find  the  magnitudes  of  all  the  stresses  and  dis¬ 
placements  occurring  during  the  bending  of  a  three-layer  axlsymmetrlcal 
shell,  regardless  of  the  assumptions  E_  ■  **  or  E„  /  It  is  necessary 
to  assign  three  boundary  conditions  on  each  edge,  as  In  the  case  of  a 
homogeneous  axlsymmetrlcal  shell: 

a)  for  a  rigidly  embedded  edge 

(84) 


or 

tc-u— P— 0; 

b)  for  a  shifting  pinched  edge 

a>— {$— Ni“0; 

c)  for  an  edge  supported  on  hinges 

u—Afj—iO 

or 

d)  for  a  freely  supported  edge 

e)  for  a  free  edge 
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Pig.  7.  For  the  determina-  Pig.  8.  Notation  for  a  con- 

tion  of  the  load  function.  ical  shell. 

dsi  ■  Rxd  0; 

dsa  «  rd  q>; 

Ra  «  s  cot  ©. 
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